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1.  Overview 

The  proposed  work  aims  to  design,  fabricate  and  test  a  new  class  of  miero-channcl  reactors, 
combining  benefits  of  (i)  micromachining  for  realizing  complex,  two-dimensional  radial 
distribution  patterns  for  heat  transfer  and  process  intensification,  and  (ii)  extrusion  for  versatile 
fabrication  of  scalable  micro-channel  networks  over  a  range  of  thermal,  mechanical  and  catalytic 
properties  not  currently  achievable  by  current  mieroreaetor  fabrication  methods.  This 
hybridization  of  two  separate  technologies  will  produce  a  new  class  chemical  reactor  capable  of 
coupling  three  or  more  separate  kinetic  and  transport  processes  within  a  single  cartridge-based 
unit,  thus  providing  a  breakthrough  in  the  design  of  portable  fuel  reformers  capable  of  efficiently 
converting  logistics  fuels  to  hydrogen.  Product  hydrogen  can  then  be  employed  by  next- 
gcncration  electrochemical  engine  (fuel  cell)  systems  at  efficiencies  significantly  greater  than 
conventional  combustion  engine  systems. 

Reforming  of  logistics  fuels  to  hydrogen  for  subsequent  utilization  in  a  fuel  cell  system  is 
comprised  of  multiple  physical  and  chemical  processes,  including  fuel  vaporization,  contaminant 
removal  and  high-tcmperature  reforming  stages.  The  need  for  autothcrmal  operation  (in  the 
absence  of  external  heat  addition,  common  to  portable  applications)  typically  requires  additional 
exothermic  reactions  (e.g.  combustion)  to  provide  sufficient  heat  to  maintain  required  reforming 
temperatures  and  supply  vaporization  heat.  The  ceramic  microchanncl  network  developed  by  the 
PI  provides  a  powerful  tool  for  coupling  all  of  these  separate  transport  and  kinetic  processes 
within  one  single  cartridge -based,  highly  efficient  reforming  unit. 

2.  Technical  Objectives 

The  proposed  research  aimed  to  investigate  the  influence  of  radial  and  axial  heat  transport  in 
micro-channel  networks,  for  developing  a  new  class  of  integrated,  portable  logistics  fuel 
reformers.  The  proposed  timeline  of  research  activities  and  objectives  is  presented  in  Table  I,  and 
was  organized  to  produce  the  following  specific  deliverables,  targeted  for  the  end  of  each 
funding  period. 

•  Months  1-3:  Development  of  multi-fluid  model  and  simulation  strategics  for  process 
integration  schemes;  construction  of  experimental  apparatus.  (Period  One). 

•  Months  4-15:  Construction  and  testing  of  acrylic-packaged  prototype  systems  for  heat 
exchange  between  two  and  three  fluids,  and  comparison  with  multidimensional  modeling 
of  solid-phase  heat  conduction  (Period  Two). 

•  Months  16-27:  Construction  and  testing  of  silicon-packaged  monolith  for  heat  exchange 
between  two  and  three  non-reacting  flows,  and  comparison  with  multidimensional 
modeling  of  solid-phase  heat  conduction  (Period  Three). 

•  Months  28-36:  Demonstration  of  the  silicon-packaged  monolith  for  two-  and  three 
coupled  processes,  and  comparison  with  multidimensional  modeling  (Period  Four). 
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Table  I:  Project  Timeline 


Formulate  Multidimensional  Model,  Solution  Slralcgy 
Design  Prololype  Dislnbutors  (Acrylic) 

Phase  One  Milestone 
Construct  Prololype  Distributors  (Acrylic) 

Heal  Transfer  Experiments  in  Mini-channel  Systems 
Comparison  vsith  Multi-dimensional  models 

Phase  Two  Milestone 

Re-Design  Dislribulors,  Refine  Fabrication  Plan 
Construct  MF.MS  (Silicon)  Distributors 
Perform  High-Tcmpcralurc  Heal  Transfer  Experiment 
Comparison  wilh  Mulli-Dimcnsional  Models 

Phase  Three  Milestone 

Washcoat  Reforming  Catalysts 
Intcgralcd  Reforming  Experiments 
Comparison  wilh  Multi-Dimensional  Models 

Phase  Four  Milestone 


5.  Progress  Statement  Summary 

Funding  of  the  proposal  effort  was  initiated  07/01/2007.  As  of  10/31/2010,  all  objectives  have 
been  completed.  The  primary  outcomes  of  the  funded  effort  are: 

•  Development  of  a  new  approach  to  creating  thermally-integrated  ceramic  mieroreaetors. 

•  Experimental  demonstration  of  axial  ( 1  -D)  self-insulation  via  materials  selection. 

•  Theoretical  validation  of  axial  (1-D)  self-insulation  using  multi-fluid  models. 

•  Experimental  demonstration  of  radial  (2-D)  thermal  integration  via  architecture  selection 

•  Establishment  of  silieon-MEMS  construction  procedures  for  further  reductions  in  system 
weight. 

This  research  effort  resulted  in  the  successful  development  of  a  new  class  of  compact,  cost- 
effective  ceramic  mieroreaetors  for  thermal  integration.  The  combination  of  materials  selection 
and  integration  architecture  enables  unparalleled  manipulation  of  three-dimensional  thermal 
profiles  for  enhancing  heat  integration  in  microchcmieal  systems.  Experimental  efforts 
demonstrated  the  ability  to  realize  axial  self-insulation  with  stable  thermal  gradients  in  excess  of 
300°C,  with  radial  thermal  profiles  in  excess  of  100°C.  In  the  absence  of  external  insulation  of 
any  kind,  hydrogen  yields  of  30%  were  achieved.  Further  refinements  in  system  insulation, 
eombustion  flame  mitigation  and  distributor  designs  are  expected  to  be  capable  of  hydrogen 
yields  in  excess  of  50%  under  autothermal  (self-powered)  operation. 
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In  summary,  this  work  provides  the  groundwork  for  the  continued  development  and  refinement 
of  thermally-integrated  systems  to  achieve  further  breakthroughs  in  system  efficiency  and  energy 
densities.  Detailed  presentation  of  all  research  results  is  provided  in  Section  5. 

4.  Research  Outcomes  &  Discussion 

4.1  Objectives  for  Period  One  (06/01/07  -  09/30/07): 

•  Formulate  multi-fluid  model,  solution  strategy 
.otype  distributors  (Acrylic) 

;eve  mment 

Multi-fluid  mouv.is  were  developed  for  analyzing  the  influence  of  fluid  and  solid-phase 
thermal  properties,  reaction  kinetics,  thermodynamics,  and  packaging  conditions 
(adiabatic/insulat;"g  or  isothcrmal/uninsulating)  upon  thermal  efficiency  and  reaction 
nroercssio’-  d  endo-  and  exothermic  process  streams.  In  this  section,  model  derivations 

are  detailed. 

A  general  pseudo-homogeneous,  onc-dimcnsional  steady-state  model  is  developed  to 
describe  plug-flow  in  two  parallel  microchanncls  separated  by  an  impermeable  wall,  shown  in 
Figure  1  .The  model  assumes  incompressible  flow  and  neglects  mass  transfer  limitations  between 
catalyst  and  fluid  and  internal  diffusion  limitations  within  the  catalyst  phase.  A  unique  fluid  is 
supplied  to  each  channel  at  fixed  velocity  (u^i.i)  and  inlet  temperature  (T  1-1,2)  in  cither  co-  or 
counter-current  configuration.  Each  fluid  experiences  chemical  reaction  and/or  heat  transfer  with 
the  divief  wall  The  dividing  wall  in  turn  exchanges  heat  with  both  fluids,  and  also  with  its 
surround  via  packaging  conditions  at  either  end  of  the  device. 


muz//. 
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Figure  1:  Schematic  of  simplified  one-dimensional  model  describing  heat  transfer  between 
two  parallel  microchannels  with  solid-phase  heat  conduction  in  axial  direction 

This  preliminary  study  concentrates  solely  upon  the  effect  of  solid-phase  heat  dispersion  on 
heat  transfer  efficiency;  the  further  assumption  of  negligible  axial  dispersion  is  also  employed. 
The  steady-state  energy  balance  governing  each  fluid  thus  accounts  for  heat  generated  via 
chemical  reaction,  heat  transfer  between  fluid  and  the  solid  wall.  (Equation  1.1) 


3/54 


B.A.  Wilhite 


Final  Report,  2007  YFA 


10/31/2010 


ha. 


pmFp, 


{L  -T)- 


Kjc: 

-exp 

'-Ea,' 

PM.Cp, 

RT, 

(1.1) 


The  mass  balance  for  fluid  i,  accounting  for  a  simplified  irreversible,  power  law  kinetic 
reaction  in  the  absence  of  mass  transfer  limitations  and  with  axial  convection  in  either  co-  or 
counter-current  operation  is  given  by  Equation  1.2: 


dC 

±—±  =  rA  -~kiUC,  exp 
dz 


-Ea. 


RT 


(1.2) 


The  resulting  one-dimensional  model  describing  the  dividing  wall  accounts  for  axial  thermal 
conduction  and  heat  addition/removal  by  the  two  fluid-phases  (Equation  3): 


(13) 


The  influence  of  packaging  upon  heat  dissipation  is  modeled  by  appropriate  selection  of 
boundary  conditions  for  the  solid-phase.  This  analysis  considers  the  two  ideal  eases  of  (l) 
isothermal  and  (li)  adiabatic  packaging  of  the  microchanncl  network.  Isothermal  boundary' 
conditions  assume  that  the  wall  temperature  is  fixed  at  either  extrema.  For  counter-current 
systems,  wall  temperatures  arc  assumed  equal  to  the  corresponding  inlet  fluid  temperature  at 
each  extrema.  For  co-current  systems,  wall  temperatures  arc  assumed  equal  to  the  weighted- 
average  temperature  of  the  two  fluids  at  each  extrema.  Adiabatic  conditions  assume  that  the  first 
derivative  of  the  wall  temperature  at  each  extrema  is  fixed  and  equal  to  zero.  In  general,  the 
isothermal  and  adiabatic  boundary  conditions  may  be  written  as: 


Tw  =  fixed  @  z~0y  L 


(1.4a) 


— -  «  0  @  z  =  0,  L  (1.4b) 

dz 


Temperature,  concentration  and  position  within  the  microchanncl  network  arc  rendered 
dimensionless  as  follows: 


T,(z)-Tmm  .y  c,(z). 

u  =  — - - — -  — — , 


T  -T 

max  mm 


(1.5) 


The  dimensionless  temperature  0,  is  defined  in  terms  of  the  maximum  and  minimum 
achievable  temperatures  within  the  heat-exchanger  system.  The  resulting  dimensionless 
expressions  for  M  fluids  exchanging  heat  through  a  single,  radially  uniform  wall  arc 


dY_ 

dz 


-Day"  exp 


-Ea, 

RT; 


(IT’) 
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From  Equation  1.2’,  the  kinetics  of  the  ith  reacting  fluid  is  described  by  Damkohler  number, 
(Dai),  dimensionless  activation  coefficient  (a,)  and  dimensionless  heat  of  reaction  (y,)-  Fluid- 
solid  heat  transfer  is  described  by  the  number  of  transfer  units  (NTU,),  while  the  solid-phase 
axial  conduction  is  described  by  the  conduction  parameter  (CP,),  following  previous  analysis  of 
non-reacting  fluids  by  Peterson  (1999)  and  Sticf  and  coworkers  (1999).  Value  ranges  for  each 
parameter  were  estimated  based  upon  typical  microchannel  dimensions  and  different  solid-phase 
material  properties.  Fluid-solid  heat  transfer  coefficients  are  estimated  assuming  gas  phase,  while 
fluid-solid  heat  transfer  coefficients  are  estimated  from  the  correlations  of  Sieder  and  Tate 
(Gcankoplis,  2003).  Details  of  parameter  space  estimation  arc  presented  in  Table  1.1. 


Isothermal  boundary  conditions  allow  pinning  the  wall  temperature  at  either  boundary,  i.c. 
the  packaging-microehanncl  interface,  allowing  free  heat  exchange  between  the  microchannel 
solid-phase  and  the  packaging  layers.  These  conditions  correspond  to  the  limiting  case  of 
packaging  with  high  thermal  conductivity  materials  (e.g.  silicon)  which  allow  thermal 
equilibration  between  the  packaging  layer  and  the  inlet  fluid.  Solution  requires  determination  of 
heat  flux  between  packaging  layer  and  solid-phase  at  the  lower  boundary  or  initial  condition 
(z=0)  that  satisfies  target  isothermal  temperature  at  the  upper  boundary  (z=L).  For  the  case  of 
eounter-current  operation,  this  upper-bound  wall  temperature  is  equal  to  the  inlet  temperature  of 
counter-current  flow.  For  co-current  operation,  the  upper  boundary  wall  temperature  is  equal  to 
the  mass  averaged  temperature  of  the  two  exiting  fluids.  Adiabatic  boundary  conditions  assume 
zero  heat  flux  between  the  packaging  and  the  microchannel  network,  in  turn  allowing  solid-phase 
wall  temperatures  to  be  driven  purely  by  heat  exchange  between  the  two  fluids.  These  conditions 
correspond  to  the  limiting  case  of  packaging  with  thermally  insulating  materials.  Solution 
requires  determination  of  lower  boundary  wall  temperature  to  satisfy  the  zero  flux  condition  at 
the  upper  boundary,  independent  of  fluid  temperatures. 

4.1. 1.1.  Modeling  Analysis,  Case  1:  Microscale  Heat  Exchanger 

The  simplest  ease  of  the  general  model  is  the  steady  state  heat  exchange  between  two  fluids 
in  absence  of  any  chemical  reaction  (a„  Dai  =  0).  The  resulting  system  of  equations  is  cast  solely 
in  terms  of  the  number  of  transfer  units  (NTU)  and  fluid  conduction  parameter  (CP). 

^-NTUfa-e)  (1.6) 

^■•NTU2(ew-Q2)  (1.7) 
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d% 

dz 2 


CP 


(1.8) 


The  heat  exchanger  effectiveness,  e,  is  calculated  from  the  resulting  temperature  profile  as 
follows: 


=  02 


(1.9) 


Additionally,  for  the  isothermal  ease,  the  ratio  of  conductive  (solid-phase)  to  convective 
(fluid-phase)  heat  losses  can  be  obtained  from  Equation  1.10,  as  originally  shown  by  Peterson 
(1999): 


k  =  CP 


(dejdzl 


(1.10) 


Results  obtained  for  the  ease  of  isothermal  packaging  conditions  are  presented  in  Figure  2 


iai  16) 


Figure  2:  Heat  Exchanger  Simulation  Results,  for  Isothermal  Packaging  Conditions:  (a) 
effectiveness  and  (b)  conduction  losses  as  a  function  of  NTU  and  CP;  corresponding 
thermal  profiles  for  NTU=5  at  (c)  CP=0.1,  (d)  CP=1,  (e)  CP=10 

Heat  exchanger  efficiencies  and  conduction  losses  over  a  range  of  NTU  are  presented  for  the 
eases  of  low,  intermediate  and  high  solid-phase  thermal  conductivities  in  Figure  2a,  with 
corresponding  temperature  profiles  of  the  solid  and  both  fluid-phases  in  Figure  2-e  for  NTU=5. 
As  seen  in  Figure  2a,  the  effectiveness  is  a  weak  function  of  the  conduction  parameter 
However,  as  CP  increases,  substantial  increases  in  conduction  heat  losses  occur,  For  this 
isothermal  case,  the  effectiveness  is  held  at  high  values  despite  increasing  conduction  losses  with 
increasing  CP,  because  the  wall  temperature  is  pinned  by  external  conditions.  These  results 
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clearly  demonstrate  the  benefit  of  employing  intermediate  thermal  conductivity  wall  materials  to 
maintain  rapid  heat  transfer  while  minimizing  eonduetive  heat  losses  via  packaging. 

Results  obtained  for  the  ease  of  adiabatie  packaging  conditions  are  presented  in  Figure  3. 
Efficiencies  arc  presented  for  the  cases  of  low,  intermediate  and  high  solid-phase  thermal 
conductivities  in  Figure  3a.  Heat  exehanger  effectiveness  is  significantly  redueed  with  increasing 
CP,  despite  the  removal  of  eonduetive  heat  losses  (via  adiabatic  packaging  conditions).  This  ean 
be  understood  by  examining  the  corresponding  temperature  profiles  of  the  solid  and  fluid-phases 
presented  in  Figure  3b-d  for  NTU=5.  As  CP  increases  the  system  approaches  the  limit  of  thermal 
equilibration  between  the  two  fluids  and  the  wall,  equivalent  to  eo-eurrent  operation,  and  in  turn 
imposing  a  maximum  effieieney  of  50%.  Conversely,  low  thermal  conductivity  materials  allow 
maintenance  of  substantial  thermal  gradients  along  the  solid-phase  axial  length  by  preventing 
thermal  equilibration. 


Figure  3:  Results  for  Adiabatic  Microscale  Heat  Exchanger:  (a)  effectiveness  as  a  function 
of  NTU  and  CP;  corresponding  thermal  profiles  for  NTU=5  at  (b)CP=0.1,  (c)CP=l, 

(d)CP=10 

In  both  eases,  the  effectiveness  is  a  weak  function  of  CP  at  low  NTU  where  the  solution  is 
dominated  by  the  NTU,  and  under  isothermal  conditions  the  effectiveness  is  dominated  by  CP  at 
large  NTU  owing  to  aeeumulatcd  conduction  losses. 

4.1. 1.2.  Modeling  Analysis,  Case  II:  Microscale  Regenerative  Combustor 

The  analysis  is  extended  to  the  ease  of  one  exothermic  reacting  fluid  exchanging  heat  with  a 
seeond  non-reaeting  fluid  via  the  solid-wall,  for  the  ease  of  an  irreversible  first  order  reaction. 
Both  fluids  are  assumed  to  enter  the  system  at  identical  temperatures  in  order  to  isolate  the 
exothermic  reaction  as  the  sole  heat  souree  in  the  system.  With  this  condition,  the  minimum 
temperature  achievable  in  the  system  beeomes  the  coolant/reaetant  inlet  temperature 
(Tmjn=T| jniet)  while  the  maximum  theoretical  fluid  temperature  is  obtained  solely  from  the 
adiabatie  temperature  rise  of  the  exothermic  reaction.  Equations  1  ’-3’  are  rewritten  as  follows: 
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dz 


=  - Da ,Fj  exp 


1  +  y,0, 


(1.11) 


^  =  NTU,(en.  -  0, )  +  DalFj  exp 

-j?-  =  NTU2(tin  -02) 

020„  NTU 


-a, 


1  +  YA 


dz  CP, 


(1.12) 

(1.13) 

(1  14) 


Equations  11-14  allow  to  study  the  effects  of  design  (NTU,  CP)  and  kinetic  parameters  (a,  y. 
Da)  alongside  both  adiabatic  and  isothermal  conditions  on  system  performance.  This  is  cast  in 
terms  of  reaction  conversion  and  system  thermal  efficiency,  the  latter  defines  as: 


n:.,~T2 


2  \jnlei 


T  -T 


A7\ 


0, 


(1.15) 


In  order  to  investigate  the  effect  of  NTU  and  CP  on  system  thermal  efficiency  and  reactor 
conversion  for  both  the  isothermal  and  the  adiabatic  packaging,  a  set  of  “mild”  kinetic 
parameters  are  initially  employed  (y=2,  ci=  10). 


Figure  4:  Results  of  Regenerative  Combustion  Simulations;  (a)  isothermal  packaging 
conditions;  (b)  adiabatic  packaging  conditions. 

Figure  4  represents  a  parametric  study  varying  CP  and  NTU  for  both  packaging  conditions. 
CP  is  varied  from  0.1  to  10  over  a  range  of  NTU’s,  to  simulate  a  wide  selection  of  solid-phase 
thermal  conductivities.  A  sufficient  value  of  the  Damkohler  number  is  employed  to  ensure 
complete  conversion  in  the  absence  of  steady-state  multiplicity  (Moreno,  et  al.  2008).  Results  in 
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Figure  4a  clearly  show  that  the  effectiveness  is  a  strong  function  of  the  conduction  parameter  for 
the  isothermal  case,  unlike  the  heat  exchanger  ease  discussed  in  Figure  2.  This  difference  can  be 
understood  by  examination  of  the  wall  boundary  conditions.  In  the  heat  exchanger  analysis,  wall 
temperatures  arc  fixed  equal  to  fluid  inlets;  i.e.  the  driving  force  for  heat  transfer  is  defined 
solely  by  the  boundary  conditions,  not  by  the  system.  For  reacting  fluid  analysis,  wall  boundary 
temperatures  arc  fixed  corresponding  to  identical  inlet  temperatures,  such  that  the  exothermic 
reaction  (occurring  within  the  system)  is  the  sole  driving  force  for  heat  transfer.  Increasing  NTU 
over  the  range  of  ~0-5  results  in  increasing  effectiveness  regardless  of  CP,  while  further 
increases  in  NTU  result  in  reduced  effectiveness,  as  conduction  losses  accumulate. 

Figure  4b  shows  the  results  for  the  effect  of  CP  and  NTU  on  the  adiabatic  heat  exchanger 
reactor.  At  low  NTU,  the  effectiveness  is  a  weak  function  of  the  conduction  parameter  in  the 
absence  of  conduction  losses  via  packaging.  However,  at  high  NTU,  CP  again  becomes  an 
important  parameter,  reflecting  the  potential  for  localized  heat  accumulation  in  the  solid-phase 
(i.e.  hot-spot  formation).  For  the  remainder  of  the  analysis  presented,  a  value  of  NTU=5  is 
assumed,  corresponding  to  a  balance  between  providing  sufficient  contacting  time  for  heat 
transfer  between  fluids  and  minimizing  conduction  losses  and/or  reaetor  size.  The  influence  of 
solid-phase  wall  conduction  upon  hot-spot  formation,  steady  state  multiplicity  and  thermal 
efficiency  at  fixed  intermediate  values  of  NTU  is  studied  in  detailed  (see  Moreno  et  al)  for  the 
cases  of  isothermal  and  adiabatic  packaging  conditions,  respectively.  In  both  eases,  sufficient 
low  solid-phase  thermal  conduction  supports  the  formation  of  a  localized  hot-spot  along  the 
solid-phase  axial  profile,  in  addition  to  the  expected  fluid-phase  hot-spot.  The  effect  of  thermal 
conductivity  on  the  system  performance  and  the  presence  of  steady-state  multiplicity  are  also 
studied.  Ignition/extinetion  behavior  is  expected,  with  important  implications  for  microsystem 
operation.  For  the  ease  of  isothermal  packaging,  the  rate  of  heat  removal  from  the  system  via 
packaging  losses  increases  (independent  of  convective  heat  removal),  while  serves  to  stabilize 
the  system  at  the  expense  of  substantial  losses.  For  the  case  of  adiabatic  packaging,  only  the 
unstable  solutions  change  with  conduction  parameter,  while  the  range  of  Da  corresponding  to 
multiple  steady-states  is  unaffected  by  CP.  This  analysis  suggests  a  design  trade-off  between 
reaetor  stability  and  competitive  thermal  efficiencies. 

It  is  important  to  note  that  for  the  case  of  isothermal  packaging  and  sufficient  CP  and  NTU 
corresponding  to  isothermal  coolant  and  solid-phase  conditions,  the  system  of  equations  reduces 
to  a  single  plug-flow  reacting  fluid  described  by  an  initial-value  ordinary  differential  equation 
which  must  yield  a  unique  solution.  The  introduction  of  boundary- value  equations  describing  the 
solid-phase  temperature  as  well  as  any  counter-current  flows  in  turn  allows  the  possibility  of 
multiple  steady-states.  That  the  present  model,  which  neglects  several  established  sources  of 
steady-state  multiplicity  demonstrates  multiplicity  due  solely  to  packaging  and  materials 
selection  clearly  shows  the  importance  of  these  two  design  parameters  in  the  design  of  thermally- 
integrated  microsystems. 

4.1. 1.3.  Modeling  Analysis,  Case  III:  Heat  Exchanger  Reactor 

Lastly,  analysis  is  extended  to  the  coupling  of  an  exothermic  reacting  fluid  with  an 
endothermic  reacting  fluid.  The  general  model  Equations  1  ’-3’  reduce  to  the  following: 
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dY, 

dz 


-Da,  F,  exp 


i  +  yA 


(1.16) 


— -  =  ±Da-X^  exp 
dz 


-ct2 

.1  +  yA 


(1.17) 


dO, 


^-/V7't/,(01,-0,)  +  Da,F1exp 


</0, 


-a, 


1  +  yA 


— f  =  NTUJew  -e2)±  Da2K  exp 

c/e 


-a. 


1  + 


</20...  ATD, 


K-et)^(e..-e2) 


dz 2  CP, 


CP, 


(1.18) 

(1.19) 

(1.20) 


Thermal  efficiencies  are  calculated  based  on  the  heat  released  from  the  exothermic  reaction 
that  is  utilized  to  (i)  drive  the  endothermic  reaction  and  (ii)  as  sensible  heat;  as  follows: 

£  =02  +  X2y2  (1.21) 

Results  from  Figure  5  show'  that  system  efficiency  and  stability  are  complex  functions  of  the 
effect  of  conduction  parameter.  Under  isothermal  conditions,  increasing  CP  results  in  losses  of 
thermal  efficiency  (at  sufficient  values  of  the  Damkohler  number  corresponding  to  complete 
conversion  of  the  exothermic  reacting  fluid)  due  to  conduction  losses  via  packaging,  as  observed 
in  the  microscale  regenerative  combustor  case.  Under  adiabatic  conditions,  the  solutions  of  very 
low  and  very  high  values  of  Damkohler  number  are  independent  of  the  conduction  parameter, 
and  thermal  efficiencies  are  consistently  higher  owing  to  the  absence  of  conductive  losses  via 
packaging.  In  both  cases  increasing  CP  results  in  multiple  steady-states  for  both  reacting  fluids. 
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(a)  (b) 


Figure  5:  Simulation  Results  for  Co-Current  Heat  Exchanger  Reactor,  (a,c)  isothermal 

packaging;  (b,d)  adiabatic  packaging. 

Similar  analysis  for  the  case  of  counter-current  operation  is  presented  in  Figure  6. 
Comparison  of  results  in  Figure  5  (for  co-current)  and  Figure  6  (for  counter-current)  illustrates 
the  additional  influence  of  flow  configuration  on  both  thermal  efficiency  and  multiplicity.  While 
maximum  efficiencies  are  similar  for  both  co-  and  counter-current  flow,  Damkohlcr  numbers 
necessary  for  complete  conversion  of  the  exothermic  fluid  reduce  mildly  for  the  case  of  counter- 
current  flow,  mostly  noticeably  at  low  values  of  the  conduction  parameter.  Multiplicity  behavior 
is  similar  for  both  adiabatic  cases  while  counter-current  operation  removes  steady-state 
multiplicity  in  the  isothermal  case. 


(a)  (b) 


Figure  6:  Simulation  Results  for  Counter-Current  Heat  Exchanger  Reactor,  (a,c) 
isothermal  packaging;  (b,d)  adiabatic  packaging. 
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4.1.2.  Summary  of  Modeling  Analysis 

Modeling  analysis  indicates  that  under  expected  operating  conditions  and  materials  selection, 
heat-exchanger  efficiencies  of  60  -  80%  are  predicted  for  the  silicon-packaged  ceramic  mini¬ 
channel  system.  Simulations  of  regenerative-combustion  configuration  (exothermic  reacting 
fluid,  heat  transfer  fluid)  predict  thermal  efficiencies  approaching  30%  for  isothermal 
packaging  conditions ,  and  efficiencies  approaching  65%  for  adiabatic  packaging .  These 
results  indicate  the  need  for  minimizing  heat  transfer  from  the  monolith  to  the  packaging 
materials,  either  by  minimizing  interfaeial  contact  area  or  the  use  of  insulating  gasket  materials 
during  packaging.  Simulations  of  heat-exchanger  reactor  configuration  (exothermic  reacting 
fluid,  endothermic  reacting  fluid)  predict  thermal  efficiencies  approaching  60-70?4t  albeit  with 
significant  light-off  or  ignition  behavior .  Thermal  efficiencies  are  improved,  at  the  expense  of 
magnified  ignition  behavior,  by  minimizing  heat  transfer  between  the  eordicrite  mini-channel 
network  and  the  silicon  distributors.  Overall,  simulations  predict  promising  thermal  efficiencies 
while  identifying  the  factors  influencing  steady-state  multiplicity  and  ignition/extinction 
behavior .  Simulations  justify  the  selection  of  eordicrite  ceramics  as  sufficiently  low  thermal 
conductivity  materials  for  the  proposed  effort. 
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4.2.  Objectives  for  Period  Two  (10/01/07  -  09/30/08): 

•  Construct  Prototype  Distributors  (Acrylic) 

•  Distribution  Quality  and  Heat  Transfer  Experiments  in  Mini-channel  Systems 

•  Comparison  with  Multi-dimensional  models 

4.2.1.  Construction  of  Prototype  Distributors 

The  first  prototype  distributor  has  been  completed  and  fully  analyzed  for  fluid  distribution 
and  heat-transfer  performance.  Based  upon  findings  on  distnbution  quality,  additional  pressure- 
drop  constrictions  were  incorporated  into  all  three  designs. 


Figure  7:  Checkerboard  acrylic  prototype  for  two-fluid  integration 


4.2.2.  Distribution  Quality  Measurements  Using  Mini-channel  System 
4.2.2. 1.  Residence  Time  Distribution.  Liquid-Phase  Flow 

Residence  time  distributions  were  obtained  using  the  experimental  apparatus  and  prototype 
distributor  (Figure  7b).  Residence  time  distributions  were  determined  experimentally  by  tracer 
injections  of  aqueous  blue  dye  solutions  into  the  reactor.  Absorbance  measurements  in  the 
effluent  stream  are  taken  as  a  function  of  time  using  a  UV-Vis  spectrophotometer  and  flow¬ 
through  cuvette  at  a  constant  wavelength  of  430  nm.  The  absorbance  vs.  time  profiles  were 
obtained  for:  i)  Flow  through  a  single  cordicrite  channel,  li)  5x5  microchanncl  reactor  with  tracer 
injection  through  12  channels,  iii)  5x5  microchanncl  reactor  with  tracer  injection  through  13 
channels. 
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Figure  8:  Schematic  representation  of  the  RTD  analysis 


The  Bccr-Lambert  law  states  that  the  absorbance  of  a  solution  is  directly  proportional  to  the 
solution's  concentration.  Thus  UV/V1S  absorbance  readings  can  be  used  to  determine  the 
concentration  of  a  solution.  Under  the  assumption  that  The  Bccr-Lambert  law  applies,  the 
absorbance  profiles  can  be  directly  used  to  calculate  the  Residence  Time  Distribution  function, 
E(t),  (sec  Equation  2.1)  as  follows: 


F(t)m — - 

u  />)*  /»* 


(2.1) 


Where  A  is  the  measured  absorbance,  C(t),  is  the  concentration  of  the  blue  dye  solution  with 
time  and  F,  is  the  constant  of  proportionality  and  is  calculated  based  on  the  molar  absorptivity 
and  the  length  of  light  path  through  the  cuvette.  This  factor  F  is  assumed  to  be  constant  through 
the  experiment. 

The  first  and  second  moments  of  the  RTD,  or  the  mean  value  and  the  variance,  respectively 
are  calculated  as  follows: 


-fitE(t)dt 


(2.2) 


^ (2-3) 

These  two  parameters  arc  used  to  characterize  the  RTD  and  to  calculate  the  dispersion 
coefficient.  Assuming  Danckwcrts  boundary  conditions  for  an  open-open  system,  the  Peclct,  Pe, 
(Bodestein,  Bo)  number,  defined  as  the  ratio  of  the  rate  of  transport  by  convection  to  rate  of 
transport  by  diffusion  or  dispersion,  Pe  =  Ul/D,  can  be  calculated  from  the  following  correlation 


Pe 


(2A) 
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Table  2:  Results  of  RTD  analysis,  with  calculated  Peclet  (Dispersion)  values 


tm(s) 

Variance 

(s2) 

Pe 

1/Pe 

Flow 

rate 

(ml/min) 

Fluid 

velocity 

(m/s) 

Effective 

Diffusivity 

(m2/s) 

Single  channel 

8.4 

40.348 

5.881 

0.17 

8 

0.021 

0.00054 

5x5co-12 

23.219 

507.157 

4.167 

0.24 

2.5 

0.00667 

0000238 

5x5co-13 

17.3377 

426.259 

3.183 

0.314 

2.31 

0.006154 

0.000288 

Liquid  RTD 


Time  ($) 

Figure  9.  Liquid  RTD  Response  for  single-channel  reference  case  and  for  12-  and  13- 
channels  in  cocurrent  feed  configuration  in  5x5  distributor 

From  Figure  9,  it  ean  be  seen  that  there  is  some  dispersion  attributed  to  the  5x5  distribution 
prototype  being  tested  that  causes  the  dye  solution  pulse  to  broaden  as  it  moves  down  the 
channels.  The  Absorbanee-time  curve  for  the  5x5  distributor  has  a  long  tail  and  this  causes  a 
large  spread,  therefore  large  variance  that  contributes  to  very  different  Pe  numbers  (see  Table  2) 
compared  to  the  RTD  through  a  single  ehannel.  As  shown  in  Error!  Reference  source  not 
found.,  values  for  1  /Pe  for  all  three  experiments  are  in  the  range  of  0.1  -  1,  indicating  that  mild 
dispersive  effeets  should  be  expeeted  in  the  experimental  system.  In  all  three  eases,  RTD 
responses  follow  theoretical  shapes  eommiserate  with  diffusive/dispersive  mixing,  and  not  from 
flow  maldistribution  (which  is  expeeted  to  produee  multiple  peaks,  each  corresponding  to  a 
unique  flow  path).  Thus,  the  present  analysis  indicates  that  the  distributor  prototypes  are 
capable  of  uniform  distribution ,  with  mild  contributions  to  dispersive  mixing  due  to  additional 
flow  lengths . 

4. 2. 2.1.  Residence  Time  Distribution.  Gas-Phase  Flow 
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Similar  procedure  as  above  was  performed  to  calculate  the  Residence  Time  Distribution 
Function  for  the  gas  phase.  The  residence  time  distribution  was  determined  experimentally  by 
gas  injections  of  Argon  into  the  reactor  under  steady-state  Flelium  flow.  Argon  partial  pressure 
was  measured  in  the  effluent  stream  as  a  function  of  time  using  a  Mass  Spectrometer.  The 
pressure  vs.  time  profiles  were  obtained  for:  i)  Flow  through  a  single  cordieritc  channel,  ii)  5x5 
microchannel  reactor  with  gas  injection  through  12  channels,  iii)  5x5  microehannel  reactor  with 
gas  injection  through  1 3  channels. 
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Figure  10:  Gas  RTD  Response  for  single-channel  reference  case  and  for  12-  and  13- 
channels  in  co-current  feed  configuration  in  5x5  distributor 

Additional  measurements  were  performed  in  order  to  compare  the  RTD  for  a  5x5 
microehannel  network  assembled  with  a  top  plate  designed  with  channel  diameter  of  1/16”  and  a 
plate  designed  with  channel  diameter  of  1/32”  for  additional  pressure  drop  constriction. 


Gas  RTD 
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Figure  11:  Gas  RTD  Response  for  a  5x5  acrylic  distributor:  a)  top  plate  channel 
diameter=l/16”,  b)  additional  pressure  drop  constriction,  top  plate  channel  diameter=l/32’ 
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Table  3:  Results  of  gas  residence-time  distributor  analysis,  including  calculated  Peclet 

(Dispersion)  values 


tm  (S) 

Variance 

(s2) 

Pe 

1/Pe 

Flow  rate 
(ml/min) 

Fluid 

velocity 

(m/s) 

Effective 

Diffusivity 

(m2/s) 

Single  channel 

36.3 

103.721 

28.936 

0.0345 

10 

0.0267 

0.000137 

5x5col2-l/16” 

38.17 

421.09 

9.757 

0.102 

10 

0.0267 

0.0273 

5x5col3-l/16” 

35.27 

213.33 

14.81 

0.0675 

10 

0.0267 

0.0180 

5x5col  2-1/32” 

36.52 

335.2 

10.884 

0.0919 

10 

0.0267 

0.0245 

5x5col3-l/32” 

32.77 

270.44 

10.867 

0.092 

10 

0.0267 

0.0245 

From  the  RTD  results  presented  in  Figure  10  and  Figure  1 1  we  eould  verify  for  all  the  eases 
studied,  that  there  is  some  dispersion  along  the  path  of  the  gas,  which  may  be  attributed  to 
laminarity  of  flow  and  diffusion.  Sinee  the  RTD  eurve  shows  only  one  peak  it  can  be  eoneluded 
that  the  distributor  does  not  present  channeling  and/or  gas  crossover  within  channels.  The 
difference  in  areas  and/or  effective  diffusivities  ean  be  attributed  to  the  dilution  faetor  used  on 
the  experiments.  For  the  single  channel  reactor,  1  ml  of  Argon  was  injeeted  in  a  stream  of  He/N: 
at  a  flow  rate  of  10  ml/min.  For  the  5x5  networks,  1  ml  of  Argon  was  injeeted  in  a  stream  of 
He/N:  at  a  flow  rate  of  1 20  ml/min. 

Additional  experiments  were  performed  after  incorporating  pressure  drop  constrictions  at  the 
distributor-monolith  interface;  these  constrictions  are  1/32"  ID  (as  compared  to  the  original  1/16" 
outlet  ID),  thus  providing  an  identical  limiting  pressure-drop  region  for  each  flow  path  This,  in 
turn,  is  expected  to  provide  greater  flow  uniformity  aeross  the  distributor.  Comparing  the  results 
obtained  the  5x5  prototype  with  and  without  the  pressure  drop  eonstrietion,  it  ean  be  verified  that 
the  top  plate  designed  with  smaller  diameter  (1/32")  provides  a  better  flow  distribution  Error! 
Reference  source  not  found,  presents  a  summary  of  ealeulated  dispersion  coefficients  obtained 
for  gas-phase  flows. 


4.2.3.  Heat  Transfer  Experiments  in  Mini-channel  Systems 

Initially  the  3x3  aerylie  prototype  was  used  to  verify'  model  results  for  Case  1  (heat  transfer 
between  two  non-reaeting  fluids).  Cold  water/hot  water  system  was  used  to  measure  the  heat 
transfer  effectiveness  (e)  and  to  monitor  the  temperature  profile  along  the  axial  direction  Overall 
heat  transfer  effectiveness  under  non-reaeting  conditions  was  ealeulated  by  monitoring  inlet  and 
outlet  temperatures  for  both  hot  and  eold  flows.  Axial  thermal  profiles  for  eounter-eurrent  flow 
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are  obtained  by  embedding  multiple  thermocouples  within  the  ceramic  monolith  wall.  Data 
acquisition  is  automated  through  a  PC  equipped  with  Labview  software. 

In  order  to  verify  the  model  results  for  Case  1,  multiple  runs  of  different  water  flow  rates 
were  tested  for  counter-current  configuration.  The  heat  transfer  effectiveness  was  calculated 
from  Equation  7,  based  on  the  temperature  data  collected  from  inlet  and  outlet  flows.  From 
Figure  12,  it  can  be  seen  that  heat  transfer  effectiveness  of  about  20-40%  were  achieved  for  both 
cases.  Increasing  the  water  flow  rate  for  both  configurations  consistently  decreases  the 
effectiveness.  For  instance,  the  effectiveness  calculated  at  a  flow  rate  of  8  ml/min  corresponding 
to  NTU=2.1  is  of  about  35%,  increasing  the  flow  rate  to  15  ml/min  for  a  NTU=1.4  decreases  the 
effectiveness  to  25%.  Also,  it  is  noticeable  that  the  configuration  “hot5cold4”  (hot  water  flowing 
through  5  channels  and  cold  water  through  4  channels)  is  more  thermally  efficient  than 
“hot4cold5”  configuration. 
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Figure  12:  Heat  transfer  effectiveness  for  counter-current  3x3  microchannel  network  for 
a)hot  water  flow  through  5  channels  and  b)  hot  water  flow  through  4  channels 
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Figure  13:  Axial  temperature  profiles  for  a)  Cold  and  b)  Hot  fluid  flowing  in  5-channel 

network  (3x3  prototype). 

Heat  transfer  experiments  for  cold/hot  water  system  in  the  5x5  checkerboard  acrylic 
prototype  were  also  performed  to  validate  the  model  results  of  case  1:  Microscale  Heat 
Exchanger.  Effectiveness  results  from  Figure  14  how  that  increasing  the  flow  rate  decreases  the 
effectiveness  as  shown  for  the  3x3  prototype  in  Figure  1 2  (from  50%  to  43%).  Moreover,  the  5x5 
acrylic  distributor  prototype  is  certainly  more  efficient  than  the  3x3  network  (effectiveness  -0.4- 
0.5  for  5x5  prototype  vs  -0.2-0. 3  for  3x3  prototype).  This  can  be  explained  by  the  extra  layer  of 
channels  that  the  5x5  networks  holds  that  adds  more  contact  area,  therefore  more  efficient  heat 
transfer. 

The  addition  of  pressure  drop  constrictions  by  reducing  the  diameter  of  the  distributor  top 
plate  from  1/16"  to  1/32  not  only  docs  enhance  the  flow  distribution  as  discussed  on  the  previous 
section,  it  also  improves  the  heat  transfer  in  the  device.  Effectiveness  results  from  Figure  14 
confirms  the  last  statement  by  showing  an  increase  of  about  10%  with  the  decrease  in  channel 
diameter. 
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Figure  14:  Effectiveness  calculation  for  3x3  and  5x5  acrylic  prototypes 


4.2.3.  Comparison  with  Multi-Dimensional  Models 

In  the  modeling  analysis,  Case  1:  Microscale  Heat  Exchanger,  the  most  important 
dimensionless  parameters  arc  the  Number  of  Transfer  Units  (NTU)  and  the  Conduction 
Parameter  (CP),  as  defined  by  Equations  2. 5-2. 6.  The  parameter  values  are  calculated  from  the 
experimental  conditions  for  the  hot-eold  water  system  and  using  the  appropriate  thermal 
properties  of  both  w  ater  and  ceramic  wall,  as  detailed  in  Tabic  4. 

halxt  kw 

NTU  =  -L-^r~ ,  CP  = - (2. 5-2. 6) 

P,Cp<  P.u.Cp, 
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Table  4:Expected  parameter  values  for  Modeling  Case  1:  Microscale  Heat  Exchanger 


Design  properties 
diameter  (m)  2.50E-03 

length  (m)  0.15 

surfaee  area  (m:)  1.50E-03 

area  6.25  E-06 

volume  (m3)  9.375E-07 

a  (area/volume)  1600  00 

Thermal  properties 
p(kg/m3)  1000 

Cp  (J/kgK)  4187 

p(Pa.s)  0.0010 

k(w/mk)  0.58 

kw  (w/mK)  2 

Prandlt(20°C)  7.01 

flow  rate(ml/min) 

5 

10 

15 

velocity  (m/s) 

0.0133 

0.0267 

0.04 

*((  s) 

11.25 

5.63 

1 

3.75 

Reynolds 

33.33 

66.67 

100 

h(w/m2K) 

678.9 

855.38 

979.2 

NTU 

2.91 

1.83 

1.4 

CP 

3. 58e-5 

1. 79e-5 

I.19e-5 

As  shown  in  Table  4,  the  experimental  values  for  the  NTU  arc  within  the  range  studied  in  the 


model  simulations.  However,  the  experimental  CP’s  are  lower  compared  to  the  simulated  values. 
To  be  able  to  present  a  straight  comparison  between  experiments  and  model  prediction,  a  low 
limit  value  for  CP  was  used  to  obtain  model  predictions  for  both  adiabatie  and  isothermal 
packaging  conditions.  This  is  due  to  highly  stiff  ODE  system  that  is  unstable  under  very  low 
eonduetion  parameter  values. 
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F’igure  15:Hot/cold  water  system,  flow  rate  =  10  ml/min.  Model  calculations  for  NTU=1.83 

and  CP<0.004 

Simulated  and  experimental  axial  temperature  profiles  for  the  “cold  flow"  stream  are 
presented  in  Figure  15.  The  results  clearly  show  a  substantial  agreement  between  simulated 
model  results  for  both  adiabatic  and  isothermal  packaging  conditions  and  the  experimental  data 
for  Case  1:  Microscale  Heat  Exchanger.  The  model  predicts  £=0.476  and  £=0.498  for  the  ideal 
eases  of  adiabatic  and  isothermal  conditions  respectively  compared  to  the  experimental  value  of 
e=0.4. 


4.5.2.4.  Summary 

Comparison  of  modeling  results  with  experimental  data  for  the  case  of  heat  transfer  between 
two  non-reacting  fluids  was  performed  using  aqueous  solutions  in  low-temperature  acrylic 
prototypes.  Experimental  results  agreed  with  model  predictions  of  thermal  efficiencies  ranging 
from  40  -  50%.  Comparison  of  fluid  thermal  profiles  with  model  predictions  indicates  good 
agreement  between  theory  and  the  prototype  system. 
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4.3.  Objectives  for  Period  Three  (10/01/08-05/31/10): 


•  Re-dcsign  distributors,  refine  fabrication  plan 

•  Washcoat  reforming  catalyst 

•  Integrated  reforming  experiments 

•  Construct  mems  (silicon  distributors) 

4.3.  Brass-Based  Distributors  for  Microchannels 

4.3.1.  Construction  of  Intermediate-Prototype  Brass  Distributor 

Three  intermediate  brass  distributor  prototypes  were  produced  for  the  purpose  of  verifying 
high-temperature  packaging  methods,  catalyst  introduction  and  obtaining  preliminary  data  ahead 
of  Si-MEMS  construction  of  final  distributors.  The  distributor  designs  were  based  upon  the 
original  nine-channel  checkerboard  two-fluid  distributor  discussed  in  the  original  project 
proposal. 

4.3. 1.1.  Investigation  of  High-Temperature  Ceramic  Sealants 

Multiple  high-temperature  ceramic  sealants  were  investigated  for  irreversible  packaging  of 
the  twenty  five  -channel  cordientc  mini-channel  network  to  the  brass  distributors.  A  graded 
adhesive  joint  was  used  to  bond  the  cordieritc  monolith  to  the  brass  plates  due  to  gross  difference 
in  Coefficient  of  Thermal  Expansion  (CTE)  (Brass  CTE  ~19xl0'6  1/°C,  cordierite's  CTE  -0.5- 
0.6x1 0'6  1/°C)).  First  each  substrate  was  coated  with  the  adhesive  that  best  matched  its  CTE.  A 
layer  of  Aremco-Bond  598-A  (CTE  11.7x10 6  1/°C)  was  applied  to  the  brass  plates;  likewise  the 
cordierite  monolith  was  coated  with  Ceramabond  618N  (CTE  0.59x1  O'6  1/°C).  Then,  a  third 
adhesive  with  an  intermediate  CTE  (Pyro  duel  598-A  CTE  1 1 .7x10 6  1/°C)  was  used  to  bond  the 
parts  together.  After  numerous  trials,  it  was  discovered  that  the  adhesives  chosen  were  porous 
and  did  not  have  the  consistency  required  for  the  brass  plate  and  the  ceramic  monolith  to 
successfully  bond.  The  final  mixture  was  grainy  with  large  particulates,  thus  inappropriate  for 
providing  a  precision  seal  within  the  tight  machine  tolerances  of  the  distributor  plate. 

After  this  first  attempt,  ceramic,  fiber,  alumina  and  metallic  based-type  adhesives  from 
Cotronics  were  evaluated.  901  Fiber  based  Alumina ,  989  Alumina  Ceramic ,  950  Durabond 
Metallic  Adhesive,  907  GF  Resbond  Adhesive  and  Sealing  Putty.  The  lack  of  success  with  any 
of  the  adhesives  was  attributed  to  the  lack  of  surface  contact  between  the  monolith  and  the  brass 
plate,  material  cracking  after  exposing  the  bonded  material  to  high  temperatures,  glue 
consistency  degraded/weaken  when  in  contact  with  water  and  ultimately  the  monolith  would 
break  away  from  the  brass  plate.  Based  upon  these  results,  distributor  design  was  altered  to 
enable  the  use  of  reversible  compression-gasket  seals. 

4.3. 1.2.  Investigation  of  Graphite  Compression  Gaskets 

After  exploring  a  wide  variety  of  high  temperature  sealants  for  irreversible  packaging  of  the 
cordieritc  network  to  the  brass  distributors  and  experiencing  the  issues  explained  above,  a 
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modified  brass  distributor  relying  upon  a  flat,  machined  graphite  compression  gasket  with 
thickness  of  1/16”  for  realizing  leak-free  and  reversible  sealing  of  mini-channel  network  to 
distributors  was  built. 

The  mieroreaetor  is  comprised  of  two  brass  distributors,  sealed  to  either  end  of  the  extruded 
ceramic  mieroehannel  network  via  compression  chuck.  Distributors  were  assembled  from 
individually  machined  brass  plates,  and  subsequently  laminated  to  form  the  completed 
distributor,  capable  of  directly  addressing  each  mieroehannel  with  either  a  combustion 
(methanol,  oxygen)  and  reforming  (methanol,  steam)  mixture. 

For  this  study  three  intuitive  distribution  distributors  are  identified  in  order  to  investigate  the 
influence  of  radial  distribution  upon  thermal  integration  and  self-insulation  for  the  ease  of 
hydrogen  production  from  methanol  via  endothermic  steam  reforming  coupled  with  exothermic 
methanol  oxidation.  Figure  16  presents  cross-section  illustrations  of  each  integration  scheme 
employed  in  the  present  study;  the  designs  are  as  follows: 

•  Distributor  A:  checkerboard  pattern  coupling  two  separate  volumes  (combustion,  steam 
reforming)  while  maintaining  uniform  radial  temperature. 

•  Distributor  B:  annular  pattern  coupling  two  separate  volumes  (combustion,  steam 
reforming)  while  providing  radial  self-insulation. 

•  Distributor  C:  annular  pattern  coupling  three  separate  volumes  (combustion,  steam 
reforming,  insulation)  while  providing  radial  self-insulation.  For  the  present  study,  this 
third  volume  is  kept  empty  in  order  to  focus  the  study  upon  the  potential  for  additional 
radial  self-insulation  in  the  absence  of  flow-folding  (e.g.,  effluent  cool-down,  pre¬ 
heating)  benefits. 


(a)  (b)  (c) 


SR 

1 

i 

0 

t  SR  1 

l  1 

00 

E 

SR  |l  SR 

! II 

SR  ! 

_ J 

0 

1  SR 

i 

0 

SR 

El 

SR 

0 

0 

0 

nsi 

_ i 

SR 

B 

r* 

L  _J 

SR 

0 

SR 

0 

!  sr] 

SR 

S 

0 

0 

j 

_ i 

SR 

□ 

r 

! » 

1 _ 

0 

!  sr 

0 

SR 

Q] 

SR 

0 

0 

0E 

u 

SR 

Ur 

1 

SR 

L_j 

0 

SR 

i 

0 

SR 

SR 

SR 

SR 

L 

001 

i 

E 

- . ;j 

l  J 

Figure  16:  Cross-section  illustrations  for  (a)  Distributor  A:  Checkerboard  pattern  for 
coupling  two  separate  volumes,  (b)  Distributor  B:  Annular  pattern  for  coupling  two 
separate  volumes  and  (c)  Distributor  C:  Annular  pattern  for  coupling  three  separate 

volumes 
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4.3. 1.3.  Verification  of  Catalytic  Washcoating  Methods 

Three  distinct  5x5  microchannel  networks  for  coupling  methanol  combustion  and  steam 
reforming  were  prepared  in  accordance  with  each  of  the  three  distributor  distributors  shown  in 
Figure  16.  For  distributor  A  (Figure  16a),  12  channels  were  individually  washcoated  with 
combustion  catalyst  and  the  remaining  13  channels  packed  with  reforming  eatalyst.  For 
distributor  B  (Figure  16b),  9  channels  were  washcoated  with  combustion  catalyst  and  the 
remaining  16  channels  packed  with  reforming  catalyst.  For  distributor  C  (Figure  16c),  a  single 
channel  was  washcoated  with  combustion  catalyst,  8  channels  were  packed  with  reforming 
catalyst,  and  the  remainder  were  left  empty.  In  all  cases,  dense  fired  cordierite  honeycomb- 
monoliths  with  a  cell  density  of  72  cells  per  square  inch  (CPS1)  were  supplied  from  Rauschert 
and  used  for  construction  of  microchannel  networks.  Preliminary  results  showed  these  substrates 
were  found  to  have  a  slight  porosity  which  resulted  in  a  subtle  flow  crossover  between  reforming 
and  eombustion  channels.  In  order  to  minimize  crossover,  cordierite  substrates  were  initially  dip- 
coated  with  a  nanopowder  alumina  layer  as  follows:  Nanopowder  alumina  slurry  was  prepared 
by  mixing  nano-sized  y-alumina  powder  (Sigma-Aldrich),  eolloidal  alumina  dispersion  (20  wt%, 
Aldrich),  distilled  water  and  methanol  (J.T.  Baker)  in  the  ratio  of  1 :0.5: 1 0.5 : 1 0.5  by  mass  (Kim 
et  al.,  2009).  Cordierite  monolith  was  dipped  into  the  slurry  and  then  ultra  sonicated  at  room 
temperature  for  1  hour.  This  procedure  was  repeated  three  times  while  allowing  sample  to  dry  at 
room  temperature  for  1  hour  between  sonication  steps.  Coated  sample  was  then  calcined  at 
650°C  for  4  hours  with  a  heating  and  cooling  rate  of  5°C/min.  This  pre-treatment  was  found  to 
remove  any  measureable  flow  crossover.  Combustion  catalyst  (1  wt%  Pt/AliOj  powder  from 
Sigma-Aldrich)  washcoating  solution  was  prepared  in  identical  fashion  to  that  detailed 
previously  (Moreno  and  Wilhite,  2010).  Blocking  gaskets  allowed  selective  coating  of 
combustion  channels  only  during  dip-coating;  monoliths  were  partially  submerged  in 
washeoating  solution  for  approximately  one  minute,  or  until  slurry  level  reaehed  the  mid-point  of 
the  microehannel  length.  Exeess  solution  was  removed  by  pressurized  gas,  and  the  resulting 
system  ealeined  at  450°C  for  4hr  (heating  rate  3°C/min  and  eooling  rate  of  3°C/min).  This 
coating  strategy  guarantees  that  half  of  the  combustion  channel  length  is  free  of  catalyst. 
Resulting  catalytic  coatings  of  thicknesses  200  -  300  pm  were  verified  by  scanning  electron 
microscopy  (SEM)  analysis  to  be  uniform  and  crack-free.  (See  Figure  17) 


Figure  17:  SEM  images  of  1  wt%  Pt/AhOj  catalyst  deposited  monolith 
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4.3. 1.4.  Verification  of  Catalytic  Packing  Methods 

For  steam  reforming,  ceramic  microchannels  were  packed  with  commercial  catalyst 
based  upon  CuO/ZnO/AliOj  formulation  (Puristar  R3 1 5  from  BASF).  Original  catalyst  is  in  the 
form  of  5  mm  pellets.  These  pellets  were  crushed  and  sieved  between  500mm  and  700mm  (25- 
35  mesh),  and  sieved  particles  packed  by  hand  into  individual  channels.  Nickel  mesh  (lOOp) 
(Alfa  Aesar  #44128)  was  inserted  into  reformer  inlets  and  outlets  to  act  as  retaining  screens. 
F.ach  channel  contained  approximately  600  mg  of  catalyst. 

4.3. 1.5.  Distributor  Construction 

Brass  plates  were  selected  for  the  construction  of  each  distributor,  owing  to  suitable 
operating  temperatures  and  case  of  machining  using  in-house  facilities.  Four  2”  x  2"  brass  plates 
were  machined  with  two-dimensional  patterns  on  cither  side,  according  to  Figure  18.  The  three 
top-most  plates  were  of  thickness  1/8”,  to  minimize  size  and  weight;  the  bottom-most  plate  was 
1”  thick  in  order  to  accommodate  conventional  1/16”  nominal  pipe  thread  (NPT)  fittings  for 
fluidic  connections  to  the  experimental  apparatus.  All  prototype  features  and  fluidic  connections 
were  1/16”,  with  the  notable  exception  of  1/32”  holes  used  on  the  top-most  plate  (at  the 
distributor-ceramic  microehanncl  interface).  These  uniform  pressure  restrictions  act  to 
compensate  for  any  variations  in  pressure  drop  associated  with  individual  flow  paths,  this 
strategy  has  been  previously  employed  in  microchemical  systems  with  complex  distribution 
networks  to  ensure  uniform  flow  distribution  (Amador  ct  al.,  2004;W'ada  ct  al. ,  2006).  Final 
distributor  assembly  was  accomplished  by  compression  scaling  the  four-plate  stack  with  silicon 
gaskets  between  each  plate.  Silicone  rubber  gaskets  (McMastcr-Carr)  were  also  used  for 
compression  scaling  cither  face  of  the  ceramic  microehanncl  network  to  brass  distributors. 

4.3. 1.6.  Microreactor  Assembly 

The  microehanncl  network  was  assembled  by  metal  plates  attached  to  each  brass  distributor 
block.  One  end  (metal  platc/brass  distributor  block)  is  fixed  and  the  other  end  is  compression 
scaled  to  extruded  ceramic  mierochannel  network  by  adjustable  screws.  Silicone  gaskets  (1/16”) 
were  placed  in  the  recessed  face  of  each  distributor,  followed  by  cither  end  of  the  ceramic 
network.  This  set-up  allows  leak-free  and  reversible  scaling  of  microehanncl  network  to 
distributors.  (See  Figure  19) 
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Architecture  A 


Architecture  B 


Architecture  C 
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Figure  18:  Schematic  of  distributor  layers  for  2-D  complex  distribution  schemes:  1. 
Bottom  plate:  Fluid  A  and  fluid  B  main  inlets;  2.  Second  plate-bottom:  fluid  A  undergoes 
first  level  bifurcation;  3.  Second  plate-top:  End  of  first  bifurcation;  4.  Third  plate:  End  of 
second  bifurcation  level;  6.  Top  plate-bottom:  Second  level  bifurcation  of  fluid  B,  7.  Top- 

top:  End  of  third  bifurcation  level 

(b) 

ff  ■PBBttj 

mmmmm 
mmmmm 
mmmmm 

mmmmm 

- 131 - 


Figure  19:  Construction  of  the  5x5  microchannel  network;  (a)  distributors,  (b)  Cross- 
section  of  microchannel  network,  (c)  side-view  of  microchannel  network, (d)  final  assembly 

with  four-point  compression  chuck. 
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4.3.2.  Microreactor  Operation 

Once  assembled,  each  of  the  three  mieroreaetors  were  investigated  over  a  range  of 
combustion  and/or  reforming  flow  rates  while  varying  the  oxygen  content  of  the  combustion 
feed.  Methanol  vapor  and  oxygen  was  supplied  to  combustion  channels  by  blending  20% 
02/5%He  in  Ar  with  ultra-high  purity  Ar  and  passing  through  an  isothermal  bubbler  containing 
99.9%  purity  methanol  (J.T.  Baker).  An  equimolar  mixture  of  distilled  water  and  methanol  was 
supplied  to  the  reforming  channels  by  passing  a  dry  5%  Ni/Ar  mixture  through  a  second  bubbler 
containing  115  mL  and  285  mL  of  methanol  and  distilled  water,  respectively.  Both  bubblers, 
thermostated  by  a  heating  tapes  attached  to  temperature  controllers  were  kept  at  20°C  such  that 
feed  molar  fraction  of  methanol  was  maintained  at  13%  and  2%  in  combustion  and  reforming 
inlet  streams,  respectively.  Gas  flow  rates  were  controlled  using  Alicat  Scientific  digital  flow 
controllers.  Dry  gasses  were  supplied  by  Airgas. 

Analysis  of  reaction  products  was  provided  on  a  dry-gas  basis  by  an  Agilent  MS3000  gas 
chromatograph,  employing  Ar  carrier  gas,  such  that  N2  and  He  served  as  internal  standards  for 
computing  effluent  flow  rates  for  reforming  and  combustion  channels,  respectively.  Two 
condensers  maintained  at  273K  and  195K  were  placed  in  series  between  gas  effluent  and  the  gas 
chromatograph  in  order  to  remove  all  moisture  and  unreacted  methanol. 


Combustion  inlet  oxygea'methanol  molar  stoichiometric  ratio,  defined  as  the  equivalency 
ratio  (ER),  was  varied  from  0.5  to  1.0  (Equation  3.1)  in  the  present  study.  Methanol  conversion 
in  both  reforming  and  combustion  channels  was  calculated  from  the  output  molar  flow  rate  of 
methane,  carbon  monoxide  and  carbon  dioxide  (Equation  3.2).  Product  yields  reported  below 
were  calculated  based  on  the  theoretical  maximum  for  each  product  (Equation  3). 
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4.3.3.  Results  and  Discussion 

4.3.3. 1  Distributor  A:  Radially  Uniform  Integration  of  Two  Volumes 

The  first  distributor  distributor  allocates  a  steam  reforming  flow  amongst  13  channels  and  a 
separate  combustion  stream  amongst  the  remaining  12  channels  to  achieve  a  thermally  uniform 
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and  balanced  integration  of  the  two  processes.  Combustion  channel  flow  rate  (Recomb)  and 
composition  (ER)  was  varied  to  identify  optimum  operating  conditions.  Combustion  equivalence 
ratio  (ER)  was  varied  from  0.5  to  1.0  for  three  different  combustion  Reynolds’  numbers  of  150, 
200  and  250  corresponding  to  total  flow  rates  of  900,  1200  and  1500  seem.  Reforming  channel 
flow  rates  were  maintained  at  a  Reynolds’  number  of  250,  corresponding  to  a  total  flow  rate  of 
1600  seem  for  all  conditions  studied.  Methanol  mole  fraction  was  maintained  at  13%  in 
combustion  channels,  and  2%  in  reforming  channels  for  all  conditions  studied. 

The  combustion  reaction  was  first  ignited  with  stoichiometric  O2/CH3OH  ratio  (ER=1). 
After  ignition,  the  steam  reforming  flow  was  started.  Once  the  flame  stabilized,  the  fuel 
composition  was  reduced  from  ER=1  to  0.5  at  intervals  of  0.1  while  maintaining  a  constant  flow 
rate.  Same  procedure  is  used  for  all  three  prototypes  at  the  flow  rates  analyzed. 

Selectively  coating  half  of  the  combustion  channel  length  with  catalyst  places  fuel  ignition  at 
the  axial  mid-point;  however,  combustion  flow  rate  and  composition  as  well  as  reforming  flow 
direction  are  found  to  have  substantial  influence  upon  final  hot-spot  location.  Co-current 
operation  resulted  in  a  stable  hot-spot  location  near  the  microehannel  outlet,  while  counter- 
current  operation  resulted  in  a  hot-spot  location  at  the  axial  mid-point;  this  is  attributed  to 
additional  combustion  heat  dispersion  provided  via  the  reforming  stream.  Thus,  all  experiments 
for  distributors  A,  B  and  C  were  performed  in  counter-current  flow. 

Hot-spot  magnitude 

Hot-spot  magnitude  as  a  function  of  ER  and  Recomb  is  presented  in  Figure  20.  At  ER  =  1 ,  hot¬ 
spot  was  located  at  the  axial  mid-point  corresponding  to  the  point  at  which  fuel  and  oxidant  first 
contact  the  combustion  catalyst.  Under  these  conditions,  the  hot-spot  has  a  symmetrical  shape 
with  a  maximum  temperature  of  560(,C  while  both  distributors  remained  at  -  50°C.  As  ER  is 
reduced,  hot-spot  position  does  not  shift;  however  the  maximum  temperature  decreases  to  a 
minimum  of  260'C  corresponding  to  an  ER  =  0.5  and  Recomb  =  150.  Further  reductions  in  ER 
resulted  in  hot-spot  extinction,  regardless  of  combustion  flow  rates. 

Extent  of  Reaction 

The  influence  of  combustion  flow  rate  (Recomb)  and  equivalence  ratio  (ER)  upon  methanol 
conversion  is  summarized  in  Figure  21a-e.  Methanol  conversion  in  the  steam  reforming  channels 
remained  constant  at  Rerer  =  250  over  the  range  of  combustion  equivalence  ratios  (ER)  studied. 
Reduction  in  methanol  supply  to  the  combustion  volume  (via  ReCOmb)  resulted  in  decreased 
methanol  conversion  in  the  reforming  channels  owing  to  reduced  system  temperature,  with  this 
effect  becoming  more  pronounced  at  ER  <  0.7.  Figure  21b  shows  the  effect  of  ER  and  ReLOmb  on 
the  performance  of  the  microehannel  combustor.  As  ER  was  reduced  from  1  to  0.5,  methanol 
conversions  in  the  combustion  channels  decreased  from  98%  to  42%.  In  addition,  methanol 
conversion  decreased  mildly  with  increasing  flow  rate,  owing  to  a  reduction  in  the  residence 
time.  Overall  methanol  conversion  (Figure  21e)  reflects  a  combination  of  these  two  trends  such 
that  total  methanol  conversion  decreases  with  eombustion  ER  regardless  of  Recomb 
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Figure  20:  Effect  of  Recamb  and  ER  upon  hot-spot  magnitude  in  distributor  A 


.  Carbon  Monoxide  Formation 

The  influence  of  both  Recomb  and  ER  upon  carbon  monoxide  formation  is  summarized  in 
Figure  21g-i.  As  the  hot-spot  magnitude  increases  through  increasing  ER  from  0.5  to  1.0,  carbon 
monoxide  yields  increase  from  10%  to  30%.  This  is  attributed  to  a  change  in  the  water-gas-shift 
equilibrium  to  the  reactant  side  as  excess  combustion  heat  accumulates  within  the  reforming 
volume,  thus  favoring  carbon  monoxide  formation.  In  the  combustion  channels,  a  maximum 
carbon  monoxide  yield  of  32%  was  observed  at  ER  =  0.8  and  ReCOmb  -  250.  Further  increases  in 
ER  reduced  carbon  monoxide  yield  as  selectivity  shifted  towards  carbon  dioxide  formation.  The 
overall  carbon  monoxide  yield  trends  are  dominated  by  those  of  the  combustion  channels,  owing 
to  the  relatively  large  amount  of  fuel  supplied  to  the  combustion  volume 

Methane  formation 

The  effect  of  both  Re  and  ER  in  the  combustion  channels  on  methane  formation  is 
summarized  in  Figure  21j-l.  In  the  reforming  channels,  the  molar  concentration  of  methane  was 
almost  always  below  detection  limits,  while  combustion  channels  achieved  <3%  methane  yield 
under  all  conditions. 
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Figure  21 :  Distributor  A:  Effect  of  ER  and  Re  in  combustion  channels  on  (a-c)  Methanol 
conversion,  (d-f)  Hydrogen  yield,  (g-i)  Carbon  monoxide  formation,  (j-1)  methane 

formation 
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Hydrogen  Production 

Figure  2 1  d-f  indicates  that  an  ER  of  0.5  -  0.6  is  necessary  to  achieve  both  high  hydrogen  and 
low  carbon  monoxide  yields.  At  ReCOmb  =  250  and  ER  =  1,  the  hot-spot  temperature  is  560°C, 
which  is  sufficient  to  achieve  hydrogen  yields  in  excess  of  90%  in  the  reforming  volume.  As  ER 
is  reduced  the  hot-spot  temperature  decreases  from  560°C  to  about  300°C  at  ER  -  0.5  and  RcC0mb 
=  250,  which  favors  the  production  of  carbon  dioxide  (Figure  21g).  Reduction  in  equivalence 
ratio  (ER)  also  favors  hydrogen  production  via  partial  oxidation  in  the  combustion  volume,  as 
shown  in  Figure  2 1  e;  a  maximum  in  combustion  volume  hydrogen  yield  of  51%  was  obtained 
for  an  ER  =  0.6  at  RcCOnib  =  250. 

The  maximum  hydrogen  yield  for  the  autothermal  reformer  reflects  a  combination  of  both 
steam  reforming  and  partial  oxidation  (combustion  volume)  performance.  As  ReCOmb  is  increased, 
the  maximum  overall  hydrogen  yield  shifts  to  lower  equivalency  ratios,  reflecting  the  importance 
of  hydrogen  yields  from  the  combustion  volume  which  consumes  a  substantial  amount  of 
methanol  fuel.  At  Recomb  =  150,  maximum  hydrogen  yield  of  28%  was  observed  at  an  ER  =  0.9; 
at  Recomb  =  200,  the  maximum  yield  corresponded  to  ER  =  0.7;  at  Recomb  =  150,  maximum 
hydrogen  yield  of  30.3%  was  observed  at  ER  =  0.6.  These  results  indicate  that  maximum 
hydrogen  yields  correspond  to  operation  of  the  heat-exchanger  microreactor  as  a  combination  of 
partial  oxidation  and  steam  reforming. 

4.3. 3. 2.  Distributor  B:  Annular  Integration  of  Two  Volumes 

The  second  distributor  allocates  a  steam  reforming  flow  amongst  16  channels  and  a  separate 
combustion  stream  amongst  the  remaining  9  channels  in  a  radial  distribution  pattern 
approximating  a  conventional  shell-and-tubc  heat  exchanger.  The  system  performance  was 
investigated  for  combustion  equivalence  ratios  (ER)  of  0.5  to  1  and  reforming  flow  rates  of  1000, 
1600  and  2000  seem,  corresponding  to  individual-channel  reforming  Reynolds’  numbers  of  125, 
200  and  250,  respectively.  Combustion  flow  rates  were  maintained  at  1125  seem,  corresponding 
to  a  Reynolds’  number  of  250.  This  value  was  selected  based  upon  best  performances  observed 
with  distributor  A. 

Hot-spot  Magnitude 

The  effect  of  combustion  ER  on  hot-spot  magnitude  for  all  three  reforming  flow  rates  is 
summarized  in  Figure  20.  The  hot-spot  was  located  at  the  axial  mid-point  of  the  reactor  and 
spanned  the  majority  of  the  mierochannel  network  with  a  maximum  external  temperature  of 
430°C.  As  ER  is  reduced,  the  hot-spot  magnitude  decreased  to  300°C  at  ER  =  0.5,  with 
distributor  temperatures  remaining  at  40  -  50°C  for  all  conditions.  Furthermore,  increasing  Rercf 
decreased  hot-spot  magnitude,  owing  to  a  combination  of  greater  heat  removal  and  utilization  by 
the  reforming  stream. 
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Extent  of  Reaction 

The  influence  of  reforming  flow  rate  and  combustion  equivalence  ratio  (ER)  upon  methanol 
conversion  is  summarized  in  Figure  22a-e.  Methanol  conversion  in  the  steam  reforming  channels 
remained  >  90%  over  the  span  of  ER  values  investigated  when  operated  at  the  lowest  reforming 
flow  rate  (Reref  =  125).  For  the  case  of  Rerei  =  200,  methanol  conversions  drop  to  60%  as  ER  is 
reduced  to  0.5;  similarly  at  Rerer  =  250  conversions  reduce  to  10%  at  ER  =  0.5.  This  data 
suggests  that  a  combination  of  insufficient  partial  oxidation/combustion  heat  and  decreased 
reforming  residence  times  significantly  limit  methanol  conversion.  This  behavior  is  similar  to 
that  observed  in  Distributor  A  when  hot-spot  magnitude  was  decreased  by  lowering  combustion 
flow  rates,  suggesting  that  thermal  effects  dominate  residence-time  effects.  Likewise,  the 
influence  of  ER  on  methanol  conversion  in  combustion  channels  is  identical  to  that  observed  for 
distributor  A,  and  can  again  be  attributed  to  a  lack  of  sufficient  oxidant  at  low  ER  values. 

Carbon  Monoxide  Formation 

The  influence  of  both  Rerer  and  ER  upon  CO  formation  is  summarized  in  Figure  22g-i.  As 
hot-spot  magnitude  increases  (via  decreasing  Reref),  both  reforming  and  combustion  carbon 
monoxide  yields  increase.  This  can  be  attributed  to  the  accumulation  of  excess  reaction  heat 
shifting  ehemieal  equilibrium  in  favor  of  earbon  monoxide  production.  Increasing  ER  from  0.5 
resulted  in  an  initial  increase  in  earbon  monoxide  formation  owing  to  this  effect;  as  ER  is 
increased  beyond  0.8,  earbon  monoxide  yields  decrease  owing  to  combustion  stoiehiometry 
favoring  complete  combustion. 

Methane  Formation 

The  results  for  methane  formation  at  the  conditions  studied  arc  shown  in  Figure  22j-l.  As  was 
the  case  for  Distributor  A,  reformate  methane  yields  were  below  the  gas  chromatograph’s 
detection  limit.  Low  (<1%)  eombustion  methane  yields  were  observed  over  the  full  range  of  Rcref 
and  ER  studied. 

Hydrogen  Production 

The  effect  of  Reref  and  ER  upon  hydrogen  yield  is  summarized  in  Figure  22d-f.  At  low  flow 
rates  (Reref  =  125),  reforming  hydrogen  yield  was  not  strongly  influenced  by  decreasing  ER, 
instead  remaining  at  -  80%.  As  noted  above,  reducing  the  combustion  ER  at  high  reforming  flow 
rates  (Rcrer  =  250)  significantly  reduced  hot-spot  magnitude  and  thus  methanol  reforming  rates 
and  hydrogen  yields.  Increasing  reforming  flow  rate  strongly  influenced  hydrogen  yields  in  the 
combustion  channels.  While  a  maximum  hydrogen  yield  of  33%  was  measured  at  low  reforming 
flows  (Reref  =  125),  less  than  10%  hydrogen  yields  were  observed  when  reforming  flow  rate  was 
doubled  (Reref  =  250).  Further  increases  in  Reref  resulted  in  reduced  methanol  conversion  in  the 
eombustion  channels,  thus  low  hydrogen  yields,  and  a  highly  unstable  behavior  prone  to 
quenching  at  low  ER  values. 


33/54 


B  A.  Wilhite 


Final  Report,  2007  YFA 


10/31/2010 


x" 


100 

80 

60 

40 

20 


N® 

ox 


V 

i— 

(N 


100 

80 

60 

40 

20 


x° 

0s* 


40 

30 


0s 


1 

0.8 

0.6 

0.4 


u  0.2 
>- 


A  A  a  £  100 

♦  *  60 
1  40 

*  X°  20 

ooooo 

O  OO'OnJ’  (N 

AO 

%  X 

e 

m- 

(C) 

V 

.5.6.7.8.91 

ER 

(d)  V®  40 

JMMt  \  30 

**  1  20 

*  x"  10 

.5.6.7.8.91 

ER 

(e)  _  30 

•  •  A  0s 

»♦*»,  '$  20 

*****!  ^  10 

.5.6.7.8.91 

ER 

(0. 

* 

.5.6.7.8.91 

ER 

#(|}  £  40 

•  ♦  *  *  8  20 

+  *  0  10 

.5.6.7.8.91  > 

ER 

so  i 

.  0.8 
0.6 
5  0.4 
0.2 

.5.6.7.8.91 

ER 

( *  •  s  40 

•  f  *.  •> 30 

■  *  ’  *  X  X  20 

f+  -  u  10 

*  .  >- 

.5.6.7.8.91 

ER 

*  't 
f* 

.5.6.7.8.91 

ER 

*  1 
>8:i 

.  *  *  —  0.4 

0  0.2 
A.  .  Sh 

.5.6.7.8.91 

ER 

(1) 

••  . 

.5.6.7.8.91 

ER 

•  Re  =125  * 

ref 

.5.6.7.8.91  r  .5.6.7.8.91 

ER  ER 

Re  =200  *  Re  =250 

ref  ref 

Figure  22:  Distributor  B  :  Effect  of  ER  and  Re  in  reforming  channels  on  (a-c)  Methanol 
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A  maximum  overall  hydrogen  yield  of  26%  was  observed  at  Reref  =  200  and  ER  =  0.8;  under 
these  conditions  both  endothermic  reforming  and  exothermic  partial  oxidation  volumes  were 
producing  hydrogen.  Combustion/partial  oxidation  hydrogen  yields  were  lower  than  those 
observed  under  similar  conditions  (ReCOmb  =  250,  ER  =  0.7)  while  carbon  monoxide  yields  were 
slightly  higher,  suggesting  that  the  annular-layered  distribution  design  did  result  in  higher 
internal  temperatures  corresponding  to  the  exothermic  combustion  channels.  However,  heat 
losses  to  ambient  from  the  outer  layer  of  steam  reforming  channels  limited  overall  hydrogen 
yields  via  necessitating  relatively  high  combustion  flow  rates. 

4.3.3.3.  Distributor  C:  Annual  core  distribution  pattern 

The  final  distributor  investigated  in  the  present  work  allocates  combustion  flow  to  a  single, 
central  channel;  steam  reforming  flow  is  allocated  in  the  8  channels  adjacent  to  the  combustion 
channel,  while  a  third  heat  utilization  /  retention  volume  is  allocated  to  the  outermost  16 
channels.  For  the  present  analysis,  this  third  volume  is  scaled  in  order  to  provide  heat  retention 
while  enabling  a  direct  comparison  between  all  three  distributors  for  the  coupling  of  steam 
reforming  and  combustion.  Findings  provide  a  basis  for  indentifying  appropriate  functions  for 
this  third  volume  in  subsequent  studies,  such  as  reactant  pre-heating  or  vaporization,  or  effluent 
cool-down  or  hydrogen  cleanup  via  water-gas-shift. 

In  reducing  the  combustion  to  reforming  volume  ratio  from  3:4  to  1 :8,  conditions  for  stable 
high-temperature  operation  changed  markedly  between  distributor  B  and  C.  Initial  experiments 
varying  the  combustion  ER  indicated  that  an  ER  of  unity  was  required  to  achieve  results  that 
were  comparable  to  previous  distributors.  For  this  reason,  performance  was  studied  at 
combustion  ER  =  1  while  varying  combustion  and  reforming  flow  rates. 

4.3.3.3.1  Influence  of  Varying  Rerer  at  fixed  Recomb  =  1000 

The  influence  of  reforming  flow  rates  upon  combustion  hot-spot  magnitude  and  reformer 
performance  was  investigated  at  a  constant  combustion  flow  rate  of  500  seem,  corresponding  to 
Retomb  =  1000.  Reformer  flow  rates  of  400,  600,  800,  1000  and  1200  seem  were  employed, 
corresponding  to  RcrCi  =  100,  150,  200,  250  and  300. 

The  effect  of  Rerer  upon  hot-spot  magnitude  is  shown  in  Figure  20.  It  is  important  to  note  that 
temperatures  reported  from  thermal  imaging  correspond  to  the  microchannel  network’s  surface 
temperature,  which  was  substantially  lower  than  those  reported  for  distributors  A  &  B;  this 
suggests  that  the  additional  heat  retention  channels  in  distributor  C  are  contributing  to  the  radial 
insulation  of  the  microchannel  network.  The  maximum  hot-spot  temperature  observed  was 
300°C,  corresponding  to  the  lowest  Rercf  investigated  (Rercr  =  100).  As  Rercf  increased  to  300, 
hot-spot  magnitude  decreased  from  300°C  to  270°C.  Additionally,  increases  in  Rcref  resulted  in 
downstream  shifting  of  the  hot-spot  location  towards  the  reforming  outlet  /  combustion  inlet. 
This  is  attributed  to  increased  combustion  heat  dispersion  by  the  reforming  stream.  Figure  23a 
shows  the  effect  of  varying  Rercf  upon  methanol  conversion.  Methanol  conversion  decreased 
from  100%  at  Rercr  =  100  to  88%  at  Rercr  =  300,  owing  to  reduced  residence  times.  Rcref  was  also 
found  to  influence  methanol  conversion  in  the  combustion  volume;  as  Rcref  was  increased  from 
100  to  300,  methanol  conversion  in  the  combustion  volume  decreased  from  66%  to  63%. 
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Figure  23:  Distributor  C:  Effect  of  Reref  on  (a)  Methanol  conversion,  (b)  Hydrogen  yield, 
(c)  carbon  monoxide  formation  and  (d)  methane  formation 


Steam  reforming  carbon  monoxide  yields  remained  relatively  steady  at  -  27%  for  all  Rerct,  as 
shown  in  Figure  23e,  while  an  ER  of  unity  in  the  combustion  channels  resulted  in  <  1%  carbon 
monoxide  yields  under  all  Reref  studied.  Thus,  overall  carbon  monoxide  yields  were  observed  to 
increase  from  3%  to  8%  with  increasing  Rercf  from  100  to  300. 

Figure  23d  shows  the  influence  of  Rercf  upon  methane  formation  Methane  yields  of  <  0  4% 
were  observed  in  the  combustion  effluent  for  all  Rercf\  corresponding  to  overall  methane  yields  of 
<  0.3%. 


The  effect  of  Rercf  upon  hydrogen  yield  is  presented  in  Figure  23b.  At  Rercf  =  100,  hydrogen 
yield  in  the  reforming  volume  is  90%;  this  decreases  to  81%  at  Reref  =  300,  owing  to  reduced 
methanol  conversion  in  response  to  lower  residence  times.  Unlike  previous  experiments  with 
distributors  A  and  B,  hydrogen  generation  only  occurred  in  the  reforming  volume;  carbon 
dioxide,  steam  and  trace  amounts  of  methane  were  the  sole  products  present  in  the  combustion 
effluent. 
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4.3.3.3.2.  Influence  of  Varying  ReCOmb  at  fixed  Rer(f  =  250 

Combustion  flow  rates  of  500,  450,  400,  350  and  300  seem,  corresponding  to  ReCOmb=1000, 
900,  800,  700  and  600,  were  employed  to  ascertain  the  extent  to  which  combustion  flow  rate 
may  be  reduced  without  altering  overall  system  performance  and  thus  maximize  hydrogen 
yields. 

The  influence  of  ReCOmb  upon  hot-spot  position  and  magnitude  is  shown  in  Figure  20.  System 
temperature  decreased  from  250°C  to  170°C  as  Recomb  decreased  from  1000  to  600.  As  noted 
above,  these  maximum  observed  temperatures  correspond  to  the  outer  surface  temperature, 
which  is  expected  to  be  lower  than  the  internal  reforming  or  combustion  temperatures.  For  all 
conditions,  axial  heat  conduction  was  sufficiently  low  that  distributor  blocks  remained  at 
temperatures  <  50°C. 

The  effect  of  ReComb  on  methanol  conversion  is  shown  in  Figure  24a.  Methanol  conversion  in 
both  reforming  and  combustion  channels  decreased  with  RcCOmb>  due  to  reductions  in  reaction 
temperatures  and  thus  reaction  rates;  further  decreases  in  ReCOmb  below  600  resulted  in  unstable 
hot-spot  leading  to  extinction. 

The  effect  of  RcComb  upon  carbon  monoxide  formation  is  shown  in  Figure  24c.  It  is 
interesting  to  note  that  while  carbon  monoxide  yields  of  -  30%  arc  observed  at  RcCOmb  <  900, 
very  low  carbon  monoxide  yields  (<  10%)  were  observed  at  Rct0mb  <  700.  This  indicates  that  the 
reduced  hot-spot  magnitude  and  width  (shown  in  Figure  20)  enabled  additional  water-gas- 
shifting  to  occur;  it  is  hypothesized  that  the  hot-spot  provided  a  catalytic  cool-down  region  in  the 
reforming  channels  that  enabled  a  rapid  shifting  of  carbon  monoxide  produced  in  the  high 
temperature  hot-spot  region.  Current  research  efforts  are  aimed  at  investigating  this  unexpected 
benefit  of  the  sharp  axial  thermal  gradients  achieved  using  this  cordicritc  microchanncl  network. 

Methane  yields  were  negligible  under  all  conditions  except  for  Retomb  =  1000,  which 
corresponded  to  9%  yield  in  the  reforming  channels  at  an  observed  hot-spot  temperature  of 
250°C 

From  Figure  24b,  hydrogen  yield  in  the  reforming  volume  was  shown  to  be  strongly 
influenced  by  Recomb,  again  owing  to  reduced  reaction  temperatures  and  therefore  rates  at 
combustion  flow  rates.  While  reformer  hydrogen  yields  were  stable  and  >  80%  for  all  Recomb  > 
700,  further  increases  in  combustion  flow  rate  resulted  in  more  fuel  combustion  and  therefore 
lower  overall  hydrogen  yields;  additionally  the  higher  reforming  temperatures  achieved  at  Retomb 
>  900  resulted  in  increased  carbon  monoxide  production  Thus,  optimal  performance  for 
distributor  C  was  observed  at  RcCOmb  =  700,  corresponding  to  an  overall  hydrogen  yield  of  32% 
and  carbon  monoxide  yield  of  4%.  Under  these  conditions,  the  external  hot-spot  temperature  was 
observed  to  be  200°C,  with  a  methanol  conversion  of  90%  and  carbon  monoxide  yield  of  10%  in 
the  reforming  volume;  data  obtained  for  identical  single-channel  reforming  flow  rates  using 
distributor  B  (Figure  22)  at  a  combustion  ER  =  0.6  corresponded  to  an  observed  hot-spot 
external  temperature  of  ~300°C.  Likewise,  at  an  identical  single-channel  reforming  flow  rate 
(Reref  =  250)  in  distributor  A,  equivalent  methanol  conversions  and  carbon  monoxide  yields  were 
observed  at  a  combustion  flow  rate  of  Ret0mb  =  250  and  ER  =  0.5,  which  resulted  in  a  hot-spot 
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external  temperature  of  ~360°C.  This  comparison  of  hot-spot  magnitude  under  identical 
reforming  flow  rates,  conversions  and  yields  strongly  suggests  that  (i)  the  annular  distribution 
scheme  introduced  in  distributor  B  contributed  to  a  reduction  in  reactor  surface  temperature  of 
-~50°C,  while  (ii)  the  additional  annular  heat  retention  layer  introduced  in  distributor  C 
contributed  to  a  further  reduction  in  reactor  surface  temperature  by  -  100°C.  This  analysis 
confirms  that  the  use  of  complex,  two-dimensional  radial  distribution  patterns  arc  capable  of 
creating  three-dimensional  thermal  gradients  which  may  be  exploited  to  enhance  hydrogen 
and/or  carbon  dioxide  yields.  However,  in  the  absence  of  external  insulation,  maximum 
hydrogen  yields  were  limited  to  -30%  for  all  three  distributors. 
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Figure  24:  Distributor  C:  Effect  of  ReCOmb  on  (a)  Methanol  conversion,  (b)  Hydrogen  yield, 
(c)  carbon  monoxide  formation  and  (d)  methane  formation 

4.3.4.  Conclusions 

Results  demonstrate  the  influence  of  distribution  patterns,  individual  process  flow  rates  and 
therefore  relative  process  enthalpies  and  contribution  of  heat  dissipation  upon  resulting  thermal 
gradients  and  reactor  performance  for  the  autothermal  production  of  hydrogen  from  methanol. 
The  three  designs  are  as  follows:  Distributor  A  represents  a  checkerboard  pattern  coupling  two 
separate  volumes  (combustion,  steam  reforming)  while  maintaining  uniform  radial  temperature 
with  a  -1 : 1  volume  ratio  of  combustion  to  reforming;  distributor  B  represents  an  annular  pattern 
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with  a  3:4  volume  ratio  of  combustion  to  reforming  while  providing  radial  self-insulation,  and 
distributor  C  couples  three  separate  volumes:  a  single  ehannel  for  eombustion  catalyst,  8 
channels  for  reforming,  while  the  outer  sixteen  channels  serving  as  a  heat  retention  volume, 
providing  radial  self-insulation. 

In  the  absence  of  any  external  insulation,  distributor  A  &  B  achieved  highest  hydrogen  yields 
when  operated  as  a  combination  of  steam  reforming  and  partial  oxidation,  owing  to  near-unity 
volume  ratios;  while  distributor  C  achieved  highest  hydrogen  yields  when  operated  as  a 
combination  of  eombustion  and  steam  reforming;  owing  to  1:8  volume  ratio  whieh  approach  the 
11:1  ratio  of  reaction  enthalpies.  Thermal  imaging  analyses  demonstrated  that  all  three 
distributors  are  capable  of  axially  isolating  hot-spots  during  operation;  this  in  turn  enabled  the 
use  of  relatively  low-temperature  reversible  gasket  sealing  and  construction  materials  for 
distributors,  while  greatly  reducing  heat  losses  to  packaging.  In  addition,  comparison  of  external 
hot-spot  magnitude  for  all  three  distributors  operating  at  identical  single-ehannel  reforming  flow 
rates  and  similar  reforming  yields  strongly  suggests  that  radial  distribution  designs  are  eapable  of 
creating  substantial  (50  -  150°C)  radial  thermal  gradients  in  the  small  (5x5)  mieroehannel 
network.  Hence  results  clearly  demonstrate  how  distribution  distributors  allow  manipulation  of 
three-dimensional  hot-spots.  Even  though  findings  demonstrate  the  value  of  axial  self-insulation 
achieving  net  hydrogen  yields  of  20-32%  for  all  three  distributors,  they  also  indicate  that 
variations  in  radial  distribution  schemes  were  not  able  to  limit  radial  heat  losses  to  surroundings. 


4.4.  Fabrication  of  Silicon  Distributors 

Using  conventional  machining  technology,  ihree  5x5  (25-channel)  prototypes  fashioned  out 
of  brass  material  were  constructed  for  the  autothermal  production  of  hydrogen  from  methanol. 
While  we  were  able  to  demonstrate  the  proof  of  concept  by  using  this  technology,  conventional 
machining  presents  restrictions  on  sealing  up  and  fabrication  of  high-preeision  microscale 
structures  with  complex  patterns  as  it  is  limited  to  sizes  down  to  1/32  in.  As  pattern  complexity 
and  size  grows  from  5x5  networks  to  larger  systems  (extrudates  are  available  on  sizes  up  to  104 
channels,  conventional  machining  becomes  unreasonable.  Si-MEMS  technology  was 
investigated  for  miniaturization  (reduced  weight  and  volume)  of  the  5x5  designs. 

This  section  details  the  design,  fabrication  and  demonstration  of  the  First  generation  of 
mierofabrieated  silicon  distributors  for  packaging  ceramic  mieroehannel  networks,  in  order  to 
reduce  overall  system  volume  and  weight.  By  using  photolithography  and  chemical  etching 
techniques  employed  by  mieroeleetromeehanieal  (MF^MS)  and  mieroehemieal  systems  the  three 
distributors  detailed  above  for  autothermal  hydrogen  production  from  methanol  are  built. 

The  fabrication  strategy  consists  of  four  silicon  wafers,  containing  a  total  of  eight  layers 
comprising  the  overall  mierofluidie  network  analogous  to  technique  described  in  Chapters  3  and 
4.  Fusion  bonding  technique  allows  high-temperature,  high-pressure  sealing  between  individual 
plates  (London  et  al.,  2001;  Mirza  and  Ayon,  1999)  to  complete  the  four-wafer  stack.  The 
resulting  distributors  will  be  compression-sealed  to  a  5x5  ceramic  mieroehannel  network  for  the 
heat  integration  of  steam  reforming  and  eombustion  of  methanol  for  hydrogen  production. 
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4.4.1.  Distributor  Designs 

Distributor  designs  utilize  silicon  mierofabrieation  methods  to  fabricate  three  complex  two- 
dimensional  distribution  patterns  across  the  mieroehannel  network.  Figure  25-27  present  3-D 
illustrations  of  each  integration  scheme  employed  in  the  present  study;  the  designs  are  as 
follows: 

•  Distributor  A:  checkerboard  pattern  (Figure  25)  coupling  two  separate  volumes 
(combustion,  steam  reforming)  while  maintaining  uniform  radial  temperature. 

•  Distributor  B:  annular  pattern  (Figure  26)  coupling  two  separate  volumes  (combustion, 
steam  reforming)  while  providing  radial  self-insulation. 

•  Distributor  C:  annular  pattern  (Figure  27)  coupling  three  separate  volumes  (combustion, 
steam  reforming,  insulation)  while  providing  radial  self-insulation.  For  the  present  study, 
this  third  volume  is  kept  empty  in  order  to  focus  the  study  upon  the  potential  for 
additional  radial  self-insulation  in  the  absence  of  flow-folding  (e.g.,  effluent  eool-down, 
pre-heating)  benefits. 


Figure  25:  3-D  schematic  of  Distributor  A 


Figure  26:  3-D  schematic  of  Distributor  B 
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Figure  27:  3-D  schematic  of  Distributor  C 
4.4.2.  Fabrication  Process 

The  fabrication  strategy  requires  machining  of  eight  unique  layers  in  four  5”  DSP  (double 
side  polished)  wafers,  400  pm  thick,  using  photolithography,  deep  reactive  ion  etching  and  wafer 
bonding  as  key  technologies.  There  are  4  dies  per  wafer  (2  for  distributor  A,  1  for  each 
distributor  B  and  C).  Eaeh  wafer  undergoes  a  sequence  of  fabrication  steps  to  form  respective 
features/ehannels  before  they  are  bonded  together  to  form  the  final  distributors.  The  planar 
photomasks  referred  to  below  are  found  in  Appendix  I  along  with  a  detailed  step-by-step  flow 
chart  of  fabrication  process  (Appendix  II). 

4.4. 2,1.  Photomasks 

The  masks  for  photolithography  were  designed  using  the  LayorrtEditor  software.  Each  mask- 
design  contained  a  set  of  alignment  marks  that  will  be  used  for  proximity  contact  (front  and 
baek)  alignment  between  masks  and  for  wafer  alignment  prior  to  bonding. 

Photomasks  eonsist  of  a  5”  polished  glass  plate,  eoated  with  a  thin  film  of  ehromium  metal. 
Subsequent  layer  of  photoresist  is  then  spun  on  the  plate,  and  the  exposure  is  made  Individual 
mask  patterns  are  drawn  into  the  photoresist  eoating  of  eaeh  photomask  plates  using  the 
GCA/Mann  3600F  Pattern  Generator  (PG).  This  tool  exposes  variably  sized  rotated  reetangles 
onto  the  mask  blank  defined  by  a  set  of  precision  controlled  moveable  aperture  blades.  Pattern 
data  is  fraetured  into  a  format  such  that  the  entire  pattern  can  be  described  by  these  rotated 
reetangles.  The  plate  is  then  developed,  exposing  the  ehromium  metal  to  be  removed. 
Photomasks  development  was  performed  using  a  Hamateeh  Mask  Plate  Processor.  Photomasks 
were  then  transferred  to  a  cleaning/stripping  tank  containing  AZ  300T  solution  (4% 
tetramethylammonium  Hydroxide,  41%  N-methylpyrrolidone,  NMP,  l-methyl-2-pyrrolidinone, 
55%  propylene  glyeol,  1,2-propane-diol)  where  they  were  dipped  for  30  minutes  in  order  to 
remove  remaining  photoresist.  The  clean  and  patterned  metal-on-glass  photomasks  were  then 
rinsed  with  DI  water  and  stored  until  needed. 
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Figure  28:  Cross-section  Schematic  of  assembled  layers  for  constructing  microfluidic 

distributors 


5.4.2.  Wafers  A,  C 

Wafers  A  and  C  each  contain  a  single  pattern,  comprising  inlet  and  through-holes 
respectively,  which  is  etched  through  the  entire  depth  of  the  substrate  (See  Figure  28).  Prior  to 
thin  film  deposition,  wafers  arc  RCA  cleaned  to  ensnre  an  ion-and-contammant  free  deposition 
surface.  RCA  cleaning  process  represents  a  standardized  chemical  cleansing  procedure 
comprised  of  a  10  minute  immersion  in  1:1:5  mixture  of  NH4OH+H2O2+H2O  at  80°C  to  strip 
organic  contaminants,  followed  by  a  30  second  immersion  in  1:50  HF:HiO  solution  to  remove 
any  native  oxide  layer,  and  finally  a  10  minute  immersion  in  1:16  HCIT^OtT^O  at  80°C  to 
remove  any  remaining  10ns  from  the  silicon  surface.  The  resulting  cleansed  silicon  substrates 
were  placed  in  a  wet  oxidation  furnace  in  order  to  grow  a  2  pm  layer  of  silicon  oxide,  whieh 
acted  as  a  hard-mask  for  subsequent  patterning.  Immediately  after  removal  from  the  wet 
oxidation  furnace,  wafers  were  placed  in  a  ehcmieal  vapor  deposition  system  and  40  nm  of 
silicon  nitride  was  deposited;  this  layer  served  as  an  additional  hard  masking  layer  for  patterning. 

Application  of  the  single  photomasks  for  wafers  A  and  C  was  accomplished  as  follows:  In 
order  to  achieve  uniform  adhesion  of  the  photoresist  to  the  substrate,  the  wafer  is  first  primed 
with  hcxamethyldisiloxane  (HMDS)  using  a  liquid  primer  solution  (20%  HMDS  in  propylene 
glycol  methyl  ether  acetate).  A  7  pm  layer  of  thick  photoresist  (Shipley  SPR220)  is  coated  onto 
the  topside  of  the  silicon  wafer  by  spinning  at  a  rate  of  3000  rpm  for  1  min.  The  photoresist  is 
then  soft  baked  at  90°C  for  1  minute  on  a  hot  plate.  The  resist  is  exposed  to  Mask  Al/Cl  by 
proximity  contact  using  the  EV620  contact  aligner  tool,  exposure  time  was  45  seconds,  wafers 
were  then  baked  at  110°C  for  2  hours,  followed  by  soaking  in  a  Tctra-mcthyl  ammonium 
hydroxide  (TMAH)  developing  solution  under  mild  agitation  until  all  UV-exposcd  photoresist 
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was  observed  to  dissolve  (~  2  minutes).  Wafers  with  fully  developed  photomasks  where  finally 
rinsed  with  D1  water  and  spun  dry.  Once  the  photomask  is  applied  and  patterned,  the  underlying 
silicon  nitride  and  silicon  oxide  hard-masks  are  developed  using  chemical  etching.  Silicon 
nitride  was  etched  via  reactive  ion  etching  (RIE)  using  CF4  etchant.  Unaxis  770  deep  reactive 
silicon  etcher  is  then  used  to  etch  400  pm  deep  patterns  of  silicon  to  make  the  through  holes.  The 
underlying  oxide/nitride  layer  provides  and  etch  stop  for  this  process.  Photoresist  is  stripped  by 
double  dipping  wafers  in  a  hot-bath  (AZ300T  solution)  and  the  silicon  nitride  is  stripped  using  a 
hot  phosphoric  acid  bath.  At  this  point,  wafers  are  ready  for  bonding.  The  underlying  oxide 
layer  is  retained  until  immediately  prior  to  wafer  bonding  in  order  to  prevent  contamination  of 
the  silicon  bonding  surfaces  by  particulates  Figure  29  shows  a  schematic  representation  of 
silicon  fabrication  steps  of  wafers  A  and  C. 

Si  wafer  (loo), 
single  crystal 
DSP,  5”  x  400  pm 

<mp—KMJLar  JJaoa 

RCA  cleaning 
Grow  oxide  (2pm) 


Deposit  40  nm 
of  S13N4 


Spin  photoresisl 
lo  mask  A) /Cl 


Wafer  C 


Wafer  A 


RIK  etch  nitride 
BOE  etch  oxide 


DRIE  etch  oxide 
400  pm 


Remove  photoresist 
Remove  nitride 


Expose,  develop 
and  post-bake 

Color  key 
Silicon  wafer 
Silicon  oxide 
Silicon  Nitride 
Photoresisl 

Figure  29:Schematic  of  microfabrication  process  for  through  holes  (wafers  A,  C) 


4.4.2.3.  W  afers  B,  D 

Wafers  B  and  D  require  a  nested  mask  in  order  to  simultaneously  form  the  front-side 
bifurcation  structure  and  the  embedded  through-holes  which  provide  fluidie  connection  to  the 
backside  bifurcation  network  in  a  single  etching  step.  This  is  accomplished  by  using  the  silicon 
oxide  layer  as  a  hard-mask  (patterned  to  masks  B1  and  Dl)  in  tandem  with  a  thick  photoresist 
layer  patterned  to  mask  B2  and  D2. 
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Fabrication  steps  (schematically  represented  in  Figure  30)  are  as  follow:  Fresh  wafers  are 
first  cleaned  using  the  RCA  process  (detailed  above)  and  subsequently  a  2  pm  silicon  oxide 
coating  is  thermally  grown  and  an  outer  40  nm  silicon  nitride  film  deposited,  in  identical 
fashioned  to  A/C  wafers. 

The  silicon  nitridc/silicon  dioxide  hard  mask  is  first  patterned  as  follows:  a  Shipley  S18B 
thin  photoresist  is  spun  on  the  wafers  (3000  rpm,  1  min),  soft  baked  at  90°C  for  1  min,  exposed 
to  photomasks  for  2  seconds,  soaked  in  TMAH  for  90  seconds  or  until  all  exposed  features  arc 
dissolved,  then  hard-baked  at  110°C  for  1  min.  The  resulting  wafers  have  bifurcation  patterns 
(Bl/Dl )  patterned  in  the  silicon  nitridc/oxide  hard  mask.  A  second  thick  photoresist  is  applied  to 
pattern  the  through  holes  (B2/D2).  Photoresist  Shipley  SPR  220,  7  pm  thick  is  spun  coated  onto 
the  wafer.  Subsequently  EV620  contact  aligner  tool  is  used  to  align  previously  patterned  Bl/Dl 
and  new  mask  B2/D2.  Figure  31(a)  shows  alignment  marks  used  for  front  contact  alignment  for 
masks  Bl/Dl  with  B2/D2.  Develop  and  post  bake.  At  this  point,  the  exposed  silicon  is  etched  to 
a  depth  of  135  pm,  forming  the  through-holes  (B2/D2)  which  will  ultimately  connect  the  front 
and  back-side  bifurcation  patterns.  The  thick  photoresist  is  then  stripped,  exposing  the  silicon 
nitnde/oxide  hard  mask.  The  exposed  silicon  (now  patterned  to  the  front  bifurcation  pattern, 
mask  Bl/Dl)  is  etched  an  additional  135  pm.  This  second  etch  then  forms  the  front  bifurcation 
pattern  to  a  depth  of  135  pm,  while  completing  the  etching  of  the  second  through-hole  pattern  to 
a  depth  of  270  pm.  Wafers  are  subsequently  stripped  and  cleansed  by  double  dipping  wafers  in  a 
hot-bath  (AZ300T  solution);  the  front  of  the  wafers  arc  then  coated  with  750  nm  of  SiO:  via 
PECVD,  which  provides  an  eteh-step  for  the  third  etehing  step. 

The  wafers  are  then  flipped  and  the  backside  is  coated  with  Shipley  SPR  220  thick  resist, 
patterned  to  D3/B3  photomasks.  Underlying  silicon  nitride  and  silicon  dioxide  were  etched  via 
reactive  ion  etching  (RIE)  using  CF4  and  CFTFyCh  etchant  respectively.  Wafcrs“are  then  DR1E 
etched  until  the  embedded  through-wafer  vias  appear  in  the  backside  features,  indicating  that 
etch  has  reached  the  SiC>2  etch  stop  Then  the  resist  is  stripped  and  the  backside  and  frontside 
silicon  nitride  is  removed  via  R1E.  Backside  alignment  mars  for  aligning  mask  B3/D3  with  the 
front  side  pattern  is  shown  in  Figure  31(b).  Both  front  and  backside  alignment  marks  shown  in 
Figure  31  are  patterned  on  cither  side  of  wafer.  By  placing  these  2-point  alignment  marks  the 
rotation  error  can  be  corrected 
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Figure  30:Schematic  cross  section  of  microfluidic  distributor  fabrication  sequence 


(b) 


a/Ai 


DX/»1 

D3/W3 


Figure  31  :Alignment  marks  for  (a)  contact  proximity,  (b)  back  side  alignment 
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4.4. 2.4.  Wafer  Bonding 

The  key  to  successful  silicon  fusion  bonding  (SFB)  (Miki  et  al. ,  2003;  Schmidt,  1998)  is  to 
maintain  smooth  and  particle  free  surfaces  throughout  processing.  The  silicon  oxide  layers 
provided  a  shielding  sealing  coat  to  protect  the  original  Si  bonding  surfaces.  After  all  wafer 
etches  were  completed  (Figure  32)  and  ready  for  bonding,  silicon  oxide  layer  is  removed  using 
buffered  oxide  etchant  (BOE)  for  30  minutes.  Wafers  were  then  immediately  rinsed  with  D1 
water  and  bonded  to  minimize  particle  deposition  on  the  wafer  surface  and  therefore  maximize 
adhesion. 


Figure  32:ActuaI  microfabricated  silicon  wafers 

Wafer  contacting  was  performed  by  using  a  Suss  SB8e  Substrate  Bonder.  The  bonding 
approach  was  as  follows:  First,  all  four  wafers  were  manually  aligned  and  contacted  to  form  a  4 
wafer  stack.  The  stack  was  then  placed  under  pressure  and  heat  (500°C,  4  bars)  in  the  bonder 
under  inert  atmosphere.  Manual  alignment  was  used  because  of  the  unavailability  of  an  optical 
alignment  chuck  at  the  CNF.  In  the  future  we  will  recommend  a  mechanical  alignment  method 
to  greatly  improve  this  step.  The  initial  bonding  step  achieves  a  reversible  adhesion  between 
wafers;  therefore  in  order  to  achieve  a  permanent  annealing'adhesion  of  the  layers  to  form  the 
final  stack,  pre-bonded  stacks  are  subsequently  placed  in  annealing  furnace  for  4  h  at  1 100°C. 

The  resulting  1.6  mm  wafer  stack  is  then  die  sawed  using  a  Disco  Abrasive  System  Model 
DAD-2H/6T  die  saw.  Ideally,  4  separate  distributors  would  be  obtained  from  each  wafer  staek. 
During  this  first  mierofabrieation  experience,  near  65%  yield  was  obtained. 


Figure  33:(a)  Wafer  stack  with  4  dies,  (b)  resulting  chip 
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4.4.3.  Experimental 

The  mieroreaetor  is  comprised  of  two  silicon  distributors,  sealed  to  either  end  of  the  extruded 
ecramie  rmeroehannel  network  via  compression  chuck  (Figure  34).  Fluidie  connections  eonsist 
of  Peek  tubes  of  1/16”  diameter  attached  to  a  square  acrylic  plate  that  serves  as  an  interface  to 
properly  plaee  peek  tubes  in  inlet  and  outlets  diameters.  Aerylie  plates  were  then  glued  to  the 
silieon  distributors  by  using  a  2-part  epoxy  glue. 

For  this  study,  three  different  distribution  patterns  (Figure  25-28)  were  employed  to  evaluate 
the  flow  distribution  using  the  silieon  prototypes.  The  eeramie  mieroehannel  networks  were 
prepared  in  the  same  fashion  explained  in  Section  3.  The  mieroehannel  network  was  assembled 
by  eompression-sealing  of  silieon  distributors  to  either  end  of  the  extruded  eeramie  mieroehannel 
network  as  follows.  Graphite  gasket  sheets  (1/16”,  MeMaster-Carr)  were  plaeed  in  the  reeessed 
faee  of  eaeh  silieon  distributor,  followed  by  either  end  of  the  eeramie  network.  Compression 
sealing  was  accomplished  by  a  four-point  compression  chuck  plaeed  around  the  assembly  and 
tightened  by  hand 


Figure  34:Distributor-paekaged  eeramie  mieroehannel  network;  (a)  unassembled  portable 
reactor  system  showing  two  distributors,  eeramie  mieroehannel  network  (5x5).  Graphite 
gaskets  and  compression  ehuek;  (b)  assembled  view. 

In  order  to  evaluate  the  functionality  of  the  silieon  distributors,  and  the  overall  assembly, 
eombustion  ehannel  flow  rate  was  started  at  0.2  SLPM  then  increased  up  to  0.5  SLPM.  Thermal 
images  were  obtained  during  steady  autothermal  operation  using  an  FLIR  Thermaeam™  infrared 
eamera.  By  operating  the  reaetor  in  the  complete  absenee  of  any  insulation,  accurate  thermal 
profiles  of  the  reactor  are  obtained.  Figure  35  shows  uncalibrated  thermal  images  taken  at  three 
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different  combustion  flow  rates.  Once  flame  was  stable,  reforming  flow  was  introduced  in  the 
system.  Thermal  images  taken  at  3  different  sets  of  combustion  and  reforming  flow  rates  are 
presented  in  Figure  36. 


F'igure  35:1  hermal  images  corresponding  to  three  different  combustion  flow  rates  (a)  0.2 

SLPM,  (b)  0.4  SLPM,  (c)  0.5  SLPM 


Figure  36:Thermal  images  corresponding  to  three  different  combustion  and  reforming  flow 
rates  (a)  Fcomb=0.5SLPM,  Fref=0.5SLPM;  (b)  Fcon,b=0.5SLPM,  Fre(=0.7SLPM,  (c) 

FComb=0.5SLPM,  Frer=1.0  SLPM 
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4.4.4.  Discussion 

From  Figure  35  and  36,  flow  non  uniformity  through  the  entire  network  is  observed  as 
combustion  flame  is  forming  only  in  the  upper  channels.  As  the  flow  rate  is  increased,  flame 
spreads  out  through  the  entire  network  increasing  the  pressure  drop.  No  flow  was  observed  at 
the  outlet  distributor,  and  a  visible  leak  appeared  between  the  ceramic  end  and  the  silicon 
distributor.  This  leak  could  be  attributed  to  flow  constrictions  such  as  graphite  gaskets,  small 
channel  diameter,  misalignment  and  glue  in  the  features  that  would  block  the  channels.  In 
addition,  sealing  the  ceramic  network  to  the  silicon  distributors  was  very  challenging  because  the 
monolith  recess  was  not  deep  enough,  leaving  limited  room  for  a  gasket,  therefore,  the  monolith 
ends  need  to  be  perfectly  polished  to  achieve  a  good  contact  with  the  silicon  wafers;  but  even 
with  a  perfect  monolith  end  sealing  represents  a  challenge.  However,  these  preliminary  findings 
have  uncovered  a  critical  number  of  refinements  that  are  expected  to  result  in  similar 
performance  to  brass  distributors. 

One  important  thing  that  can  be  highlighted  from  the  thermal  images  is  that  despite  the 
flow  maldistribution,  hot  spot  is  located  in  the  center  of  the  ceramic  monolith,  keeping  the 
distributors  at  room  temperature  thus  avoiding  future  problems  with  the  glue  method  employed 
to  attach  the  fluidic  connections  to  silicon  distributors. 

Different  chips  were  evaluated  at  both  ends  of  the  ceramic  network,  but  only  one  chip 
corresponding  to  Distributor  A  was  functional.  The  rest  of  the  chips  did  not  allow  any  flow 
through  the  channels  and  it  would  require  high  pressure  to  push  flow  through  the  silicon 
distributors. 

By  ShMEMS  the  overall  weight  and  volume  of  the  distributors  have  drastically  reduced 
Sec  Figure  37  where  both  brass  and  silicon  prototypes  are  compared.  The  initial  brass  prototype 
had  a  weight  of  690  g  while  the  silieon  distributors  with  a  weight  of  10  g  exhibit  a  weight 
reduction  of  a  factor  of  98.5% 


Precision-Machined  Brass  Plates 

#•••  v 

Lithography-Patterned  Silicon  Substates 

Figure  37:Comparison  of  overall  size  between  brass  and  silicon  distributors 
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4.4.4.  Conclusions 

Microelectromechanical  systems  (MEMS)  and  silicon  technology  are  promising  technologies 
for  miniaturizing  conventional  macro-scale  distributors  into  compact  devices.  See  Appendix  111 
for  the  projected  energy  density  using  the  silicon  packaged  ceramic  microchanncl  network 

From  the  first  micromachining  experience  building  the  silicon  distributors,  a  6%  final 
packaged  device  yield  was  obtained,  as  only  one  of  the  distributor  prototypes  was  functional. 
However,  the  main  goal  of  achieving  miniaturization  of  flow  distributors  by  reducing  weight  and 
volume  was  accomplished.  More  importantly,  this  first  effort  has  identified  the  key  refinements 
to  design  and  experimental  method  for  the  continued  development  of  this  technique.  This  low 
yield  can  be  attributed  to  notable  design  challenges  for  complex  MFMS  in  addition  to  the  lack  of 
previous  operating  experience. 

Based  on  the  experience  gained  with  the  mierofabrication  of  the  first  generation  of  silicon 
distributors,  the  follow  ing  suggestions  are  made  for  next  generation  of  silicon  distributors: 

i.  Contact  all  four  silicon  wafers  together  using  alignment  pins  for  accurate  wrafcr 
positioning. 

li.  In  order  to  improve  the  bonding  between  wafers,  it  is  recommended  to  clean  organic 
residues  from  w'afer  surfaces  with  piranha  solution  (mixture  of  sulfuric  acid  (H2SO4)  and 
hydrogen  peroxide  (H2O2)),  after  BOE  stripping  and  immediately  before  contacting. 

lii.  Design  front  wafer  with  a  deeper  monolith  recess.  This  would  allow  proper  seating  of  a 
scaling  gasket  between  the  silicon  and  ceramic  monolith.  This  could  be  achieved  by  using 
a  thicker  silicon  substrate  (c.g.  2  mm) 

iv.  Design  spigots  with  1/16“  diameter  instead  of  1/32“.  This  would  reduce  the  pressure 
drop  in  the  system  and  would  allow  flow  through  the  outlet  distributor  more  easily. 

v.  Position  inlets  and  outlets  with  a  larger  distance  in  order  to  obtain  an  optimum 
positioning  and  gluing  of  fluidic  connections;  alternately,  replace  the  hard-glue  with  a 
solder-fitting  technique  (Wilhite,  2009). 

vi.  Apply  even  pressure  when  hand  tightening  compression  chuck.  A  small  scratch  on  the 
silicon  surface  may  easily  spread  over  the  w  hole  body  under  the  influence  of  mechanical 
strain. 

Further  improvements  of  the  mierofabrication  steps  will  enable  the  construction  of  more 
functional  distributors  that  will  be  employed  for  the  demonstration  of  the  heat  integration 
between  methanol  combustion  and  steam  reforming  of  methanol  for  the  autothcrmal  production 
of  hydrogen. 
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